LOW MACH NUMBER LIMIT OF THE FULL 
NAVIER STOKES EQUATIONS 
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Abstract. The low Mach number limit for classical solutions to the 
full Navier Stokes equations is here studied. The combined effects of 
large temperature variations and thermal conduction are accounted. In 
particular we consider general initial data. The equations lead to a sin- 
gular problem, depending on a small scaling parameter, whose linearized 
is not uniformly well-posed. Yet, it is proved that the solutions exist and 
are uniformly bounded for a time interval which is independent of the 
Mach number Ma G (0, 1] , the Reynolds number Re £ [1, +oo] and the 
Peclet number Pe £ [1, +oo] . Based on uniform estimates in Sobolev 
spaces, and using a Theorem of G. Metivier and S. Schochet |31| . we 
next prove that the penalized terms converge strongly to zero. It allows 
us to rigorously justify, at least in the whole space case, the well-known 
computations given in the introduction of the P.-L. Lions' book I26| . 
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1. Introduction 

There are five key physical assumptions that dictate the nature of the 
low Mach number limit: the equations may be isentropic or non-isentropic; 
the fluid may be yiscous or inviscid; the fluid may be an efficient or a poor 
thermal conductor; the domain may be bounded or unbounded; the temper- 
ature yariations may be small or large. Yet, there are only two cases where 
the mathematical analysis of the low Mach number limit is well developed: 
first, in the isentropic regime !H1 !ni El El !2Z1 ; second, for inviscid and 
non heat-conductive fluids P !S1 !H2l • 
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Our goal is to start a rigorous analysis of the general case where the 
combined effects of large temperature variations and thermal conduction 
are accounted. As first anticipated in it yields some new problems 
concerning the nonlinear coupling of the equations. 



1.1. Setting the problem up. The full Navier-Stokes equations are: 

' dtp + d\v{pv) = 0, 
(1.1) < dt{pv) + div(/9i; v) + VP = diva, 

dt{pe) + d\w{pve) + P divv = div(A;VT) + a ■ Dv, 

where p, v = {v^ , . . . ,v'^) , P , e and T denote the fluid density, velocity, 
pressure, energy and temperature, respectively. We consider Newtonian 
gases with Lame viscosity coefficients C, and r] , so that the viscous strain 
tensor a is given by 

a := 2C,Dv + rjdivvid, 
where 2Dv = Vv + (Vw)* and is the d x d identity matrix. 

Considerable insight comes from being able to simplify the description of 
the governing equations (|1.1|) by introducing clever physical models and the 
use of judicious mathematical approximations. To reach this aim, a standard 
strategy is to introduce dimensionless numbers which determine the relative 
significance of competing physical processes taking place in moving fluids. 
Not only this allows to derive simplifled equations of motion, but also reveals 
the central feature of the phenomenon considered. 

In this paper we distinguish three dimensionless parameters: 

eG(0, 1], /UE [0,1] and kG[0, 1]. 

The first parameter e is the Mach number, namely the ratio of a character- 
istic velocity in the fiow with the sound speed in the fluid. The parameters p, 
and K are essentially the inverses of the Reynolds and Peclet numbers, they 
measure the importance of viscosity and heat-conduction. 

To rescale the equations, there are basically two approaches which are 
available. The first is to cast equations in dimensionless form by scaling 
every variable by its characteristic value j28| . The second is to consider 
one of the three changes of variables: 



(1.2) 


t - 


-> e^t, 


X - 


-> ex, 


V - 


ev. 


c- 


K, 


7] - 




k- 






t - 


-> et, 


X - 


x. 


V - 


ev. 


c- 


epC, 


V - 




k - 


~^ s i^k ^ 




t - 




X - 


x/e. 


V - 


-> ev. 


c- 


epC, 


V - 


epr]. 


k - 


-> e^Kk 



See |26( I43j for comments on the first two changes of variables. The third one 
is related to large-amplitude high-frequency solutions (these rapid variations 
are anomalous oscillations in the context of nonlinear geometric optics [H]). 
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Anyway, these two approaches both yield the same result. The full 
Navier-Stokes equations, written in a non dimensional way, are: 
' dtp + d\\{pv) = 0, 



where Q := e"pa- Dv with a ^ (namely, a = 1 for the first two changes 
of variables and a = for the last one). 

Our study is concerned with the analysis of the low Mach number limit 
for classical solutions to the full Navier Stokes equations in the non- 
isentropic general case and for general initial data. In particular, the com- 
bined effects of large temperature variations and thermal conduction are 
accounted. We are interested in the limit e — > and in proving results 
that are independent of p and k. The analysis contains two parts. We 
first prove an existence and uniform boundedness result for a time interval 
independent of the parameters e , p and k . We next study the behavior of 
the solutions when e tends to . 

Many results have been obtained in the past two decades about the justi- 
fication of the incompressible limit — which is a special case of the low Mach 
number approximation. Concerning the Euler equations (// = k = 0), the 
study began in ear lies eighties with works of Klainerman and Majda j22ll23j . 
Schochet ^36,, Isozaki .18. .19, . Ukai j^, and others. As regards the isen- 
tropic Navier-Stokes equations (/i = 1, ( and r] constants, p = p{P)), 
the mathematical analysis of the low Mach number limit has come of age 
since the pioneering works. Recent progress are presented in Danchin j8ji9|, 
Desjardins and Grenier ^Hj, Desjardins, Grenier, Lions and Masmoudi jllj . 
Hoff and Lions and Masmoudi ^7\. They are also two very interesting 
earlier results concerning the group method: Grenier |14j and Schochet [371 . 

Concerning the non-isentropic Euler equations having general initial data, 
Metivier and Schochet have recently proved a couple of theorems EOl |32] 
that supersede a number of earlier results (a part of their study is extended 
in ^ to the boundary case). In particular they have proved the existence 
of classical solutions on a time interval independent of e . The aim of this 
paper is precisely to start a rigorous analysis of the corresponding problems 
for the full Navier-Stokes equations. 

The study of the low Mach number limit is a vast subject of which we 
have barely scratched the surface here. To fill in this gap we recommend 
Desjardins and Lin J2] and Gallagher J3] for well written survey papers. 
For the reader who wishes to learn more about the physics, P.-L. Lions |26j . 
Majda |2H| and Zeytounian 43, 44 are good places to start. Detailed histor- 
ical accounts of the subject can be found in along with a broad number 
of references for further reading. In connection to the stability analysis per- 
formed below, let us point out that the research of numerical algorithms 
valid for all flow speeds is a very active field ^3 123 ESI E2] • 



(1.3) 




^ dt{pe) + div(/9we) -|- Pdivu = Kdi\{kSI'T) + eQ, 
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1.2. Uniform stability. We consider classical solutions, that is solutions 
valued in the Sobolev spaces H'^(0) with s large enough, the domain D be- 
ing either the whole space M*^ or the torus T"^ . Our main result asserts that, 
for perfect gases, the classical solutions exist and are uniformly bounded for 
a time interval independent of e, and n. We mention that the case of 
general gases involved additional difficulties (see Remarks 11.41 and II .bf) and 
will be addressed in a separate paper. 

We choose to work with the unknowns P , v = {v^, . . . jV"^) and T . In 
order to be closed, the system must to be augmented with two equations of 
state, prescribing the density p and the energy e as given functions of P 
and T . Here, we restrict ourselves to perfect gases so that 

P = RpT and e = CyT, 

where R and Cy are given positive constants. We begin by rewriting 
equations (|l.Hj) in terms of (P, f , T) . Set 7 = 1 + R/Cy ■ Performing linear 
algebra, it is found that 

( dtP + v-VP + 7Pdivt; = (7 - l)«;div(A:VT) + (7 - l)eQ, 
VP 

p{dtv + V ■ Vv) + = ^ diver. 



(1.4) 

where p = P j (RT) . The equations (|1.4() are supplemented with initial data: 



t pCvidtT + V ■ VT) +Pdwv = KdiY{kVT) + eQ. 



(1.5) P\t=o = Po, v\t=o = vo and 7[j=o = %■ 

Finally, it is assumed that C , V the coefficient of thermal conductivity 
k are C°° functions of the temperature T , satisfying 

k{T) > 0, C{T) > and r/(T) + C(T) > 0. 

Notation 1.1. Hereafter, A denotes the set of adimensioned parameters: 

A:= {a = {e,p,K)\e £ (0,1], p £ [0, 1], k G [0, 1] }. 

Theorem 1.2. Let d ^ 1 and D denotes either the whole space M.'^ or the 
torus T'^ . Consider an integer s > l + d/2 . For all positive P , T and Mq , 
there is a positive time T such that for all a = (e, /u, k) £ A and all initial 
data (Pq ,Vq,Tq) such that Pq and Tq take positive values and such that 

the Cauchy problem for (|1.4|) - ()1.5() has a unique solution (P",f",T'*) such 
that {P'' - P,v'',T'' -T) £ C°{[0,T];H'+^{B)) and such that P" and T" 
take positive values. In addition there exists a positive M , depending only 
on Mq , P and T , such that 

sup sup {e~^ ||P"(t) - P||h=(b) + + 11'^" W -2111/^.(0)1 ^ M. 
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The main obstacle is that the equations leads to a singular problem, de- 
pending on the small scaling parameter e , whose linearized is not uniformly 
well-posed in Sobolev spaces. In other words, since we consider large tem- 
perature variations, the problem is linearly (uniformly) unstable (see [SJ| for 
comments on this instability). Therefore, we cannot obtain the nonlinear 
energy estimates by differentiating the equations nor by localizing in the 
frequency space by means of Littlewood-Paley operators. In particular, the 
technical aspects are different from those present in the previous studies of 
the Cauchy problem for strong solutions of (|l.lj) .^ 

Remark 1.3. i. We will prove a more precise result which, in particular, 
exhibits some new smoothing effects for div v and VP (see Theorem 12. 7|) . 

ii. General initial data are here considered, and allow for large density and 
temperature variations. The hypothesis -Pq (a^) — P = 0(e) is the natural 
scaling to balance the acoustic components, see (|1.6j) and [9l l22 | l28l l3T] . 

iii. One technical reason why we are uniquely interested in the whole space 
R*^ or the Torus T'^ is that we will make use of the Fourier transform 
tools. A more serious obstacle is that, in the boundary case, there should be 
boundary layers to analyze 0. For the Euler equations (that is fj, = k = 0), 
however. Theorem 11.21 remains valid in the boundary case |^. 

Before leaving this paragraph, let us say some words about the difficul- 
ties involved in the proof of Theorems 11.21 After changes of variables (given 
in i i2.4j) . we are led to study a mixed hyperbolic/parabolic system of non- 
linear equations of the form: 



(1.6) 



gi{e,ep){dtp + v -Vp) + ^ div v - ■^xi(ep) div(/3(^)Ve) = Tj 

I 

g2i9,ep){dtv + v ■ Vv) + -Vp - fiB2{9,ep)v = 0, 



g3i9,£p){dt0 + v-V9)+divv-KX3iep) div{f3{9)V9) = eT 



3, 



where B2 is a second order elliptic differential operator, and Tj (i = 1,3) 
are of no consequence. 

One of the main differences between the Euler equations ^ |^ and the 
full equations is the following. When k = , it is typically easy to ob- 
tain estimates uniform in e by a simple integration by parts in which 
the large terms in 1/e cancel out (see [22||^). In sharp contrast (as ob- 
served in :28 ), when k 7^ and the initial temperature variations are large, 
the problem is more involved because the penalization operator is no longer 
skew-symmetric. Several difficulties also specifically arise for proving esti- 
mates that are independent of and k. 

Another main feature of the system H1.6() is that gi and 32 depend on 6 . 
As a consequence Vgi and Vg2 are of typical size 0(1) . The system (|1.6() 



^The latter problem has been widely studied. Starting from |29| and culminating 
in |S] |7j which investigate global strong solutions in spaces invariant by the scaling 11.21 
with /i = K = 1 , following the approach initiated by Fujita and Kato. 
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does not enter into the classical framework of singular limit jl4l \'67\ because 
of this strong coupling between the short time-scale and the fast time-scale. 
That is why we cannot derive estimates in Sobolev norm by standards meth- 
ods using differentiation of the equations. In this regard, an heuristic ar- 
gument is worth emphasis: For non-isentropic fluids, the pressure depends 
on the density as well as the temperature. Thus, as the Mach number goes 
to zero, the fact the variation in pressure converges to zero does not pre- 
clude large temperature fluctuations (provided that they are compensated 
by density fluctuations of opposite sign). Yet, even small perturbations of 
the density might result in non-small velocity variations. 

We conclude this part with a remark concerning general gases. 

Remark 1.4. For general equations of state, we are led to study systems 
having the form ()1.6|) where the coefficients xi ^ind xs depend also on 9 . 
As a consequence, the singular operator is nonlinear. Some computations 
indicate that there are no estimates uniform with respect to e . 

1.3. The low Mach number limit. We now turn to consideration of the 
limit of solutions of in as the Mach number e goes to . The 

purpose of the low Mach number approximation is to justify that the com- 
pression due to pressure variations can be neglected. This is a common 
assumption that is made when discussing the fluid dynamics of highly sub- 
sonic flows. In particular, provided the sound propagation is adiabatic, it is 
the same as saying that the flow is incompressible. On sharp contrast, this 
is no longer true if the combined effect of large temperature variations and 
heat conduction is accounted. Indeed, going back to (|1.4() we compute that, 
formally, the limit system reads 



where p = P/{RT), Cp = jCy and, in keeping with the notations of 
Theorem ll.2| P denotes the constant value of the pressure at spatial infinity. 

Theorem 1.5. Fix // G [0,1] and k £ [0,1]. Assume that (_P^,f^,T^) 
satisfy (|1.4j) and 

sup sup \\e'\P%t) - P)\\hs + \\v%t)\\jjs + \\T%t) -ZWjjs < +00, 
ee(o,i] te[o,T] 

for some fixed time T > , reference states P_, T and index s large enough. 
Suppose in addition that the initial data 7||^q — T is compactly supported. 
Then, for all s' < s, the pressure variations e^^(P^ — P) converges strongly 
to in L^{0,T;Hf^^{W^)) . Moreover, for all s' < s , {v" ,T^) converges 
strongly in L'^{f),T] Hf'^^{W^)) to a limit {v,T) satisfying the limit sys- 




(1.7) 



' 7Pdivu = (7 - l)Kdiv(A;VT), 
< p{dtv + V • Vw) -|- Vvr = p div a, 
^ pCp{dtT + vVT) = Kdiv(fcVT), 




tern (fr7|) . 
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Note that the convergence is not uniform in time for the oscihations on 
the acoustic time-scale prevent the convergence of the solutions on a small 
initial layer in time. The key to proving this convergence result is to prove 
the decay to zero of the local energy of the acoustic waves. To do so we will 
consider general systems which include (|1.6j) as a special case. In particular, 
the analysis of the general systems considered below should apply for the 
study of the low Mach number combustion as described in We mention 
that, in view of it seems possible to consider the same problem for 
exterior domains (which is interesting for aeroacoustic |44) ) . Yet, we will 
not address this question. The results proved in |H21 indicate that the 
periodic case involved important additional phenomena. 

Remark 1.6. For perfect gases, the limit constraint is linear in the fol- 
lowing sense: it is of the form divfg = with Ve := v — C^'^nkVT . In 
sharp contrast, for general gases the constraint is nonlinear. Indeed, it 
reads divv = /(P, T)Kdiv(/cVT) . As a consequence, it is not immediate 
that, in this case, the corresponding Cauchy problem for (|1.7() is well posed. 

2. Main results 

We will see in ^2.41 below that one can transform the equations (|1.4j) into 
a system of the form 

gi {6, ep) {dtp + v Vp) + ^ divv - ^Bi{9,ep)e = Ti, 
''^■'^^ ' g2{0,ep){dtv + v - Vv) +^Vp - nB2{e,ep)v = 0, 

^ gs{e,ep){dte + vve) + divv- KBs{e,ep)e = eXs, 

where the unknown (p, v,6) is a function of the variables {t,x) G M x D with 
values in M x R"^ x M . Recall that D denotes either the whole space or 
the torus T'^. Moreover, the coefficients gi, i = 1,2,3, are real- valued and 
the Bi 's are second order differential operators given by: 

Bi{e,ep) :=xi{ep)dW{m^-), 

B2{e,ep) :=X2{ep)diY{2C{e)D-)+X2{ep)V{v{e)div), 2D:=V + {V-)\ 
B3{e,ep) :=X3{ep)diY{P{e)V-). 

Finally, Tj := Xi{0,ep)F{e, ^Vv) where F G C°° is such that F{0) = 0. 
2.1. Structural assumptions. 

Assumption 2.1. To avoid confusion, denote by p) G the place 
holder of the unknown (6, ep) . 

(Al) The 5i's, i = 1,2,3, are C°° positive functions of (t?, p) G 

(A2) The coefficients /? and (" are C°° positive functions of G M , and ij 

is a C°° function of ■!? G M satisfying ri{'d) + C(i?) > . 
(A3) The coefficients Xi ^ ^ = 1;2,3, are C°° positive functions of p G M. 

Moreover, for all p G M, there holds Xiip) < Xsip) ■ 
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Assumption 2.2. Use the notation d/ = {df /dd) 6'd + {df /dp) dp. There 
exists two diffeomorphisms 3 p) ^ {S{^, p), p) £ and 3 
p) ^ {-&, gi-d, p)) G M2 such that 5(0, 0) = q{0, 0) = and 

(2.2) d5 = 53d^?-5idp and dg = -^g^d^ + gidp. 

X3 

Brief discussion of the hypotheses. 

(i) The main hypothesis in Assumption 12.11 is the inequahty Xi < X3 • 
It plays a crucial role for the purpose of proving estimates (see 
the section |1] and especially the study of the simplified system (|4.5|) '). 
Moreover, given the assumption /?(??) > 0, it ensures that the operator 
Bi(9,0) — 33(6,0) [which appears in the last equation of the limit 
system given below in (|2.H|) ] is positive. Which means nothing but 
the fact that the limit temperature evolves according to the standard 
equation of heat diffusion! 

(ii) The identities given in (|2.2|) are compatibility conditions between the 
penalization operator and the viscous perturbation. The reason for 
introducing S [resp. g] will be clear in ^6.51 [resp. i jT.lj . 

2.2. Uniform stability result. Given ^ t ^ T, a normed space X and 
a function 99 defined on [0, T] x D with values in X , we denote by (p{t) 
the function D B x ^ ip{t, x) G X . The usual Sobolev spaces are de- 
noted i/'^(D) . Recall that, when D = R"^ , they are equipped with the 
norms 

where u is the Fourier transform of u and (0 := (1 + ICP)^''^ • As concerns 
the case B = T'^ , replace the integrals in ^ G M'^ by sums on /c G Z*^ . 
Let us introduce a bit of notation we use continually in the sequel. 

Notation 2.3. Given cj G M and g^ 0, set ||?i||j:/<T := ||n||j|^CT-i + H'^^Hj^ct • 

RecaU that A := {(e,/i,K)|e G (0,1], /i G [0,1], k G [0,1]}. The norm 
that we will control is the following. 



Definition 2.4. Let T > , a ^ , a = {e, n, k) ^ A and set v := + k . 

The space TC^iT) consists of these functions {p,v,9) defined on [0, T] x B 
with values in M. x x M. and such that \\{p,v,9)\\'^,t^j'-^ < +00, where 

WIT) '■— sup 

{\Mt),vit))\\H?+^ + \\m\H;;+^} + 

t6[0,T] 



H + K ||V6'||^^+i + Hi \\divv\\'jja + {iJ, + k) \\Vp\\%a dt 

Similarly, the space TLao consists of these functions {p,v,6) defined on 
with values in M x M'^ x M and such that 11 (p, t;, ^) 11.^^ < +00, where 

||(p,'(;,6')||.;^<. := ||(p,'(;)||^a+i + ||6'||^<.+i with u := ^J^^ + K. 



Remark 2.5. One may not expect at first to have to take these norms 
into account, but they come up on their own in ^ Moreover, not only the 
norms, but also the spaces depend on a (since we allow fj, = or k = 0). 

Notation 2.6. Given a normed space X and a non-negative M , we denote 
by B(X;M) the ball of center and radius M in X . 

Here is our main result. In the context of Assumptions I?!n42.2l we prove 
that the solutions of (|2.1j) exist and are uniformly bounded for a time interval 
which is independent of a G A. 

Theorem 2.7. Suppose that system (|2.1|) satisfies Assumptions \2.1\ - \2.^ . 

and let d ^ 1 . For all integer s > 1 + d/2 and for all positive Mq , there 
exists a positive T and a positive M such that for all a G A and all initial 
data in B (Til^ q] Mq) , the Cauchy problem for (|2.1|) has a unique classical 
solution in b{hI{T);M). 

The proof of this result will occupy us till ^7.21 The crucial part is to 
obtain estimates in Sobolev norms that are independent of a G A. Notable 
technical aspects include the proof of an energy estimate for linearized equa- 
tions [see Theorem I4.l-{j and the use of new tools to localize in the frequency 
space [see Propositions \\\.2\ and 1115] . With these results in hands, we begin 
in ^ by analyzing the high frequency regime. The rest of the analysis is 
devoted to the proof of low frequencies estimates. We mention that we do 
not need specific estimates for medium frequencies. 

Remark 2.8. Up to numerous changes, a close inspection of the proof of 
Theorem 12.71 indicates that, in fact: for all M > Mq > 0, there exists 
T > such that for all o G A and all initial data in B(TC^q;Mq) the 
Cauchy problem for 1)2. 1|) has a unique classical solution in B{7i^{T); M) . 

2.3. Convergence toward the solution of the limit system. We now 

turn to considering the behavior of the solutions of (|2.1|) in R'^ as the Mach 
number e tends to zero. Fix fj, and k and consider a family of solutions 
of system (jSHJ), (p^ v^,9^) . It is assumed to be bound in C([0, T]; H'' (R'^)) 
with a large enough and T > 0. Strong compactness of 9^ is clear from 
uniform bounds for dtO^ . For the sequence {p^,v^), however, the uniform 
bounds imply only weak compactness, insufficient to insure that the limits 
satisfies the limit equations. We remedy this by proving that the penalized 
terms converge strongly to zero. 

Theorem 2.9. Suppose that the system (|2.1() satisfies Assumption \2.1\ . 
Fix fi £ [0,1] and k G [0,1], and let d ^ 1. Assume that {p^,v^,9^) 
satisfy (|2.1() and are uniformly bounded in TC^^^ ^ (T) for some fixed T > 

and s > 4 + d/2. Suppose that the initial data 9^(0) converge in 

to a function 9q decaying sufficiently rapidly at infinity in the sense that 

{x)^9o G H'^{R^) for some given 6 > 2 , where (x) := (1 + x^)-^/^ . 

Then, for all index s' < s , ^ strongly in L^(0, T; i7f^^(M'^)) and 
divv" -Xi(ep^)div(/3(r)Vr) ^ strongly in L^{0,T; H^^-^{R'^)) . 
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The proof is given in ^ It is based on a Theorem of Metivier and Scho- 
chet |31) . which we describe in Theorem 18.. SI about the decay to zero of the 
local energy for a class of wave operators with time dependent coefficients. 

We have just seen that converges to 0. The following result states 
that {v^,9^) converges toward the solution of the limit system. 

Theorem 2.10. Using the same assumptions and notations as in Theo- 
rem \2.iA . converges weakly in L'^{0,T; H'^(U.'^)) and strongly in 

L'^(0,T;Hf^^{R'^)) for all s' < s to a limit {v,e) satisfying 

g2{e, 0) {dtv + vVv)+Vtt- fiB2{e, 0)v = 0, 

g3{e,o){dte + vve)-K{x3{o)-ximdHm'^o) = Q, 

for some vr which can be chosen such that Vvr G C^{[0,T]; H'^^^( 



Given Theorem 12. 9( the proof of Theorem 12.101 follows from a close in- 
spection of the proof of Theorem 1.5 in jSJ, and so will be omitted. 

2.4. Changes of variables. To understand the role of the thermodynamics 
we rewrite equations (|1.3j) in terms of the pressure fluctuations p , velocity v 
and temperature fluctuations 6 ; where p and 9 are defined by 

(2.4) P = Pe'P, T = Te^ or p:=ilog(^^), 9:=log(^j 

where P and T are given by the statement of Theorem ll.2l 

We can convert the pressure and temperature evolution equations into 
evolution equations for the fluctuations p and 9 . Starting from (|1.4I) , it is 
accomplished most readily by logarithmic differentiation {dt^xP = ^Pdt^xP 
and dt^x'T = Tdt^x9)- By so doing it is found that {p,v,9) satisfies 



(2.5) 



P{dtp + v-Vp) + — diYv- — — — diY{kTV9) = (7 - 1)Q, 
P P 

— {dtv + V • Vv) + — Vp- IX div(2CDw) + V(?? div v) = 0, 
{dt9 + vV9)+Pdivv-K div{krV9) = eQ. 



R 

Therefore, the system ()2.1() includes ()2.5I) as a special case where 

9i ■■= -, 92 ■= 93 ■= Xi ■■= X2 = X3-= P = 

where the * indicates that the functions are evaluated at {9, ep) . Moreover, 
for i = 1,3, Ti := x*iF{0, ^/JlVv) where F{9, y/JlVv) := Q is as in (Oj) . 

We easily verify that the Assumptions I2.1H2.2I are satisfied in this case. 
Hence, Theorem 11.21 as stated in the introduction is now a consequence of 
Theorem l2.7l since P and T (given by (|2.4|) ) are obviously positive functions 
and since || • ||co([o,t];H°(d) ) ^ -^11 " llw^Cr) for some constant K independent 
of a . Similarly, Theorem 11.51 follows from Theorem 12.91 and Theorem 12.101 
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3. Localization in the frequency space 

We now develop the analysis needed to localize in the frequency space. 
The first paragraph is a review consisting of various notations which serve 
as the requested background of what follows. The core of this section is 
SS I3.2M3^ in which we prove two technical ingredients needed to localize in 
the low frequency region. This will not be used before Section El and could 
be omitted before the reader gets there. 

3.1. Preliminaries. To fix matters, in this section we work on the whole 
space , yet all the results are valid mutatis mutandis in the Torus T'^ . 
All functions are assumed to be complex valued unless otherwise specified. 
The notation d always refers to the dimension, ao always refers to a real 
number strictly greater than d/2 , and h stands for a small parameter. 

From now we use K to denote a generic constant whose value may change 
from line to line in the text. It always stand for a constant independent of h 
(whenever h can play a role). Within the proofs, we use the notation A < B 
to say that A ^ KB for such a constant K . 

We will often write the Sobolev spaces i?'^(M'^) as H'^ . Given two normed 
vector spaces Xi and X2, C{Xi,X2) denotes the space of bounded oper- 
ators from Xi to X2. We denote by |Hlxi-+X2 norm and C{X) is a 
shorthand notation for C{X, X) . 

Recall the notation 

(0 := (l + leP)'^'. 
A function q belongs to the class 5"" ( m is a given real number) if g(^) is 
a C°° function of ^ G M*^ and satisfies the differential inequalities ^ 

) fo'^ ^ G M*^ and for all multi-indices a G N"^ . Such a function 
is called a symbol. With the best constant Qa as semi-norms, 5"™ is a 
Frechet space. Given a symbol q G 5"™ , the Fourier multiplier associated 
to q is given by the operator Q acting on tempered distribution u by 
Qu := qu. For m^k G M, let A™ denotes the Fourier multiplier with 
symbol (/c^)™ . Put another way: 

(3.1) A^ = (/ - A:2A)'"/^ 

where I denotes the identity operator. 

We now introduce the first of two families of operators which are used 
in the sequel to localize in the frequency space. Let j be a function 
of ^ G ffi*^ , satisfying 

o^j^i, xo = iforieKi, = for lei ^ 2, j(e) = j(-o. 

Set jh{C) = j{h^) , for ^ /i ^ 1 and ^ G M''; so that jh is supported in 
the ball of radius 2/h about the origin. Then we define Jh as the Fourier 
multiplier with symbol jh '■ 



Jh = j{hD^). 
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Let us make a series of remarks on . The operator J/j is self-adjoint 
since jh is a real-valued function. Using that jh is even, we deduce that 
JhU is real- valued for any real-valued tempered distribution u. They are 
smoothing operators, and the family { | < /i ^ 1 } is an approximate 
identity. In particular, we will often use the simple observation that for 
all r ^ 0, there exists a non-negative constant K so that for all h G (0, 1] 
and fj G M, one has 

K 

(3.2) ||<//i||//tT^//tT+r ^ -j-^ and ||-^— JhWn'^—rH'^-'^ ^ 

One reason it is interesting to assume that j has compact support is the 
following: 

(3.3) Jh = JhJch, for all ^ c < 2"^ 

3.2. A product estimate. As alluded to above, the Friedrichs mollifiers Jh 
are interesting because they are essentially projection operators (see (|3.3() ). 
On the opposite, it is also interesting to use a family of invertible smoothing 
operators. A good candidate is the family 

{ AZ := (/ - h^Ar^^ I m ^ 0, /i G (0, 1] }. 

To begin with, we discuss the properties of these operators which are related 
to the splitting of the frequency space into two parts: the first where |^| ^ 
l/h, the second where |^| > 1/h. 

The operators depend on the scale parameter h . Heuristically, 
when m ^ 0, we expect that this /i -dependence reduces to the following 
alternative 

AAT<I and {I-Jh)AT<h'^\D,r. 
As regards the case of negative powers of , we will use that 

/i-A^-<A-™ and A^™ < /. 

These statements are made precise by the following lemma. 

Lemma 3.1. Let (mi,m-2) G be such that ^ mi ^ m2 ■ Then, for all 

h G [0, 1] and all a £ , we have 

Ci A) lift'"! A ""^2 II < T 

For all m ^ and all c G (0, 1] , there exists a positive constant K such 
that for all h G [0, 1] , and all cr G M , 

(3.5) \\j,hA^\\^^^^„^K, 

(3.6) ||(/ - J,h)^:\\h.^h^-.. ^ Kh^- 

Proof. To prove these results, we check that, on the symbol level 

o<j(c/iO(/iO~™^(2c-^r, 
Q^[i-3{chi)){hir{ir"' ^h^{c)\ 
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We are now prepared to establish the following product estimates. 

Proposition 3.2. Let do > d/2, (o-i,cj2) G and {mi, 1112) G be 
such that 

(3.7) fJi + (T2 + m-i + 7712 ^ 2cro. 

Then, there exists K depending only on d,aQ,(Ti,mi such that for all h € 
[0, 1] and Ui G iJ<^o-<Ti-m, ^ 

A — mi— 7712/ \|| ^ II A ^^Tll II II A II 

This result extends to vector valued functions. 

Remark 3.3. In words, this proposition just says that the smoothing effect 
of the operators A^*" is distributive. 

Proof. The key point is that the operators A^™ are invertible. It allows us 
to derive the desired product estimates from the corresponding results in the 
usual setting nii = m2 = . To do so Proposition 13.21 is better formulated 
as follows: there exists K such that for ah h £ [0, 1] and fi £ , 

(3.8) 

ll^^mi ma { (A"Vl) (A™'/2) } || j:^cto-<ti-c.2 ^ ^ ll/l ll/f-o--i ll/2||//<^o--2 • 

The proof of this claim is based on the decomposition of each function fi 
into two pieces: its low wave number part J^fi , and its high wave number 
part {I — Jh)fi ■ This leads to four products that are handled the same way 
introducing, for c G {0, 1} and m G M , the Fourier multipliers 

e^,^^:=h-^"^AtiJh-cI). 

That is, as is easily verified from the triangle inequality, the left-hand side 
of 1)3. 8|) is less than 



0^ci,C2 5Cl 



Hence, to prove the claim (|3.8j) it suffices now to combine three ingredients: 

\\i,cimi+C2m2 \-mi-m2 \\ ^IU,II 

II" ■"■/i ^||_H'°"0-<^l-°"2 ^ II '^ll_H'°"0-'^l-<^2-<:i™.l-c2m2 , 



f^l^2 ||//a-Q-CTi-CT2-cimi-C2m2 ^ 1 1 "^1 1 1 //''O ^ -<=1 "il 1 1 ^2 1 1 /fo-Q -<T2 - C2m2 j 

^Iblln-o-. (^e{i,2}). 



The first and last inequalities follow from Lemma l3. II In order to prove the 
second one, we first recall the classical rule of product in Sobolev spaces (see 
Theorem 8.3.1. in |EI)- For (ri,r2) G , the product maps continuously 
i?''i(R'^) X H'^^iR'^) to H''{R'^) whenever 

(3.9) 7"! + r2 ^ 0, r ^ minjri, r2} and r ^ ri + r2 — d/2, 
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with the third inequahty strict if ri or r2 or — r is equal to d/2. We next 
verify that (|3.9|) appUes with 

r := (To - (Ti - (72 - cimi - C2m2, := (Tq - <7i - for i = 1,2. 

□ 

3.3. A Friedrichs' Lemma. In this paragraph we present a result which 
complements the standard Friedrichs' Lemma. To do that we first need a 
commutator estimate, stating that the commutator of a Fourier multiplier of 
order m and the multiplication by a function is an operator of order m — 1. 

Lemma 3.4. Let ctq > d/2+1 and m £ [0, +cxd) . For any bounded subset B 
of S"" and all a £ {—cro + m, do — 1] , there exists a constant K such that 
for all symbol q £ B , all f £ i?'^o(M'^) , and all v e H''{W^) , 

(3.10) WQifu) - /Qn||^._+x ^ K \\f\\jj^, , 
where Q is the Fourier multiplier with symbol q . 

Lemma 13.41 is classical. Yet, for the convenience of the reader, we include 
a proof at the end of this paragraph. 

A word of caution: the estimate (|3.1fl|) carries over to matrix valued 
functions and symbols except for one key point. Suppose that / and q are 
matrix valued. In order that (|3.in|) be true the following condition must be 
fulfilled: q{i)f{x) = f{x)q{C) for all {x,C) G M^'^ . 

Proposition 3.5. Let ao > d/2 + 1 and m G [0, 1] . For all a in the 

interval (—(To + m, ctq — 1], there exists a constant K , such that for all 
h G (0, 1] , all f £ H^'oiRd^ and all u £ H-''''{R'^) , 

(3.11) ||A(/u)-/J;,^/||^_+, ^/i™i^||/||^.J|A^('^«+'^)n||^,. 

Remark 3.6. The thing of interest here is that the precise rate of conver- 
gence does not require much on the high wave number part of u . 



Proof. Given /u € M and g G ff^(M ) , we denote by g the multiplica- 
tion operator H~f^{W^) 3 u ^ gu £ S'(M.'^). Then, Proposition 13.51 can 
be formulated concisely in the following way. There exists a constant K , 
depending only on (Tq , m and a , such that for all h £ (0, 1] , 

||Ai-™[j,,/^]A^«+'^||^,^^, < h"^K\\f\\j,., . 

The proof of this claim makes use of the division of the frequency space into 
two pieces. The low frequencies region |^| < 1/h and the high frequencies 
region |^| > l/h. We write 

where 

C7^:=Ai— [j,,/](/-4)A^°+^ 
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and = J^-ifi- The reason for the particular choice of the constant 
will be apparent in a moment. 

Firstly, we estimate Cq . We rewrite Cq as 

(3.12) = hJ^A^-"^ [h-'^iJh - I), /] JhK"^"- 

The key point to estimate Cq is to notice that { h~^{jh — l)|0</i^l} 
is a bounded family in S"" . Indeed, if g G 5"™ and q{£,) vanishes for 
small ^, then the symbols q\{-) = X^"^q{X-) belong uniformly to 5*", 
for < A < cxD . Once this is granted, Lemma 13.41 implies that the family 
{ A^"^ [h~'^{Jh I < /i ^ 1 } is bounded in CiH") . On the other 

hand { 4A^°+'' \ < h £1} is a bounded family in CiH") , see ^U^ . 
Consequently, in light of (|3.12|) . we end up with 

<h"'\\f\\H.o- 

Our next task is to show similar estimates for . As regards , the 
fact that the operators Jh are essentially projection operators is the key 
to the proof. More precisely, we use the identity (|3.3|) written in the form 
J/i(l — J^) = [recall that = Js-i^ ]. It yields 

(3.13) = A'-^Jnf {I - J'h) A^«+^ 

It turns out that the situation is even better. Let v belongs to the Schwartz 
class S . The spectrum (support of the Fourier transform, hereafter denoted 
by spec ) of {Jhf){{I — J'^jv) is contained in 

spec (Jhf) + spec (/ - J'f,)v C B (O, 2/1^^) + B (O, lOh-^f C B (O, 3/^^^)^ 
the exterior of the ball centered at of radius 3h^^ . Which results in 

A((A/)((/- J») =0, 

that is Jh{f{{I - J'h)v)) = Jh{{I - Jh)f){{I - W) . By combining this 
identity with (|3.13l) . we are left with 

= J,((/ - J,)f)' (I - 4) A^«+^ 

To estimate HC*^ we prove a dual estimate for the operator ad- 
joint. We write (C^)* as a product of two operators: 

{c^y = {/j-»-(/ - j;;)Ar+"}{(u- A)/)' (z^-^+'^a) a'-'"}, 

where z denotes the complex conjugated of z . The problem reduces to 
establishing that 



(3.14) ||((I- A)/)'(/i'^«+V,Ai-™)||^_,^^.„ < /i- 

Indeed, since the family {h-''''-''{I - J^)A^»+'" | < /i ^ 1 } is bounded 
in C{H''0;H-'') (see (EH), the estimate (nni|) implies that 

\\{C^Y\\h-^^h-^ <h'^\\f Wh^o- 
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Which in turn imphes ||C^ < /i"* ||/||/^<to ^^^d completes the proof. 

We now have to prove (|3.14|) . Let v belongs to the Schwartz class 5. 
Since o"o > we can invoke the standard tame estimate for products, 
which leads to 



(3.15) \\{I-Jh)f\\L^ ||(/i'^o+'^AAi-™)t;||_ + 



To estimate these fours terms, we use the embedding of H'^'^ ^(M*^) into 
L°°(M'^) and the bounds given in (|3.2|) . to obtain 



W - Jh)fh^ <\\{I- A)/|Ih^o-i <h\u UH^o , 

IIU-A)/|Ih-o ^II/IIh-o, 

^ ^ II ^ II /fo'0-'"-{°'0+°'~™) ~ ^ ll^lli?— ■ 

Inserting the four estimates we just proved in ()3.15p . we conclude that 

11(1 - A)f{h'^'^''JhA'-'^v) < h"^ ||/||^.„ IIHI^-.. 

This completes the proof of (|3.14j) . □ 

Proof of Lemma To avoid trivialities (a G ) , assume that m ^ 2ao — 
1 . We establish the estimate by using the para-differential calculus of 
Bony P]. In keeping with the notations of the previous proof, denotes 
the multiplication operator u ^ fu. We denote by Tf the operator of 
para-multiplication by / . Rewrite the commutator [Q, /''] as 

[Q,Tj] + Q{f'-Tf)-{f'-Tf)Q. 
The claim then follows from the bounds 

(3.16) VaGM, l|[Q,T/]||^._^._™+i ^ci(g,cT)||/||^.o, 

(3.17) Vae (-ao,ao-l], ^C2(a)||/||^.„, 

(3.18) VaGM, IIQ||//.^i/.— <sup^|(0~"'7(0h 

where ci{-,a) : S"^ M+ takes bounded sets to bounded sets. We refer the 
reader to j2] and for the proofs of the first two inequalities (see also 
Prop. 10.2.2] and [13 Th. 9.6.4'] for the proof of (HnKl) : a detailed proof 
of 1)3. 17() is given in ^| Prop. 10.2.9]). The estimate 1)3. 18p is obvious. □ 
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4. Energy estimates for the linearized system 



Many results have been obtained in the past two decades concerning the 
symmetrization of the Navier Stokes equations (see, e.g., |^ [71 1201 1^1. 
Yet, the previous works do not include the dimensionless numbers. Here 
we prove estimates valid for all a = (e,/i, k) in A, where A is defined in 
Notation 11.11 As already written in the introduction, our result improves 
earlier works ^ I22| I31j on allowing k ^ 0. Indeed, when /t = 0, the 
penalization operator is skew-symmetric and hence the perturbation terms 
do not appear in the estimate, so that the classical proof for solutions 
to the unperturbed equations holds. In sharp contrast (as observed in ^H]^. 
when K 7^ and the initial temperature variations are large, the problem 
is more involved because the singular operator is no longer skew-symmetric. 
Several difficulties also specifically arise for the purpose of proving estimates 
that are independent of and k . In this regard we prove some additional 
smoothing effects for div v and X/p . 

We consider the following linearized equations: 



(4.1) 



giicl)){dtp + v-Vp) + ^divv-'^divi(3i{<P)Ve)= fi, 



To fix matters, the unknown (j),v,9) is a function of the variables {t,x) £ 
[0,T] X B (T is a given positive real number and ID denotes either or 
T'^ ) with values in M x R'^ x R . The coefficients (j) = (f){t, x) and v = v{t, x) 
take their values in R^ and R"^, respectively (iV is a given integer). 

Parallel to Assumption 12 . 11 we make the following hypotheses. 

Assumption 4.1. Throughout this section we require gi,g2,g?, to be C°° 
positive functions of G R^ , without recalling this assumption explicitly 
in the statements. Similarly, it is assumed that /3i , /32 , and (3^ are C°° 
functions of G R^ satisfying 

(ii > 0, /32 > 0, + /32 > 0, /33 > 0. 
Our main assumption reads 
(4.2) A < /33. 

Assumption 4.2. For our purposes, it is sufficient to prove a priori esti- 
mates. It is always assumed that the unknown U := [p, v, 6) , the coefficients 
{v, (j)) as weU as the source term / := (/i, /a, /s) are in C°([0, T]; i7°°(B)) . 

4.1. Statements of the results. We establish estimate on ||f^||'>^o(y) in 
terms of the norm ||C/(0)||^o of the data and norm of the source term /. 
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Recall from Definition 12.41 that 

nom-= sup {mv)it)\\j^, + \\m\\m} 

te[o,T] 



+ 



K ||V0||^i + ^ ||Vu||^i + K ||divt;||^2 + + k) ||Vp||^2 dt 



with u := y/^ + K. Recall the following notation we use continually in the 
sequel: given a G M and ^ ^ we set ■= |M|j:/<t-i + Q \\'\\h<^ ■ 

Theorem 4.3. There is a smooth non- decreasing function C from [0, +oo) 
to [0, +oo) such that for all a G A, all T £ (0, 1] , all coefficients v and (j) , 
and all {U,f) satisfying ()4.1() . the norm ||C^||7^0(t) satisfies the estimate 

l^o(r) ^ C(i?o)e^^(«) \\Um^o^^ + C{R) r \\{fuf2)\\Hl,, + m\Hi dt. 



where v := fi + k and 
(4.3) Ro:= Wmho^im 



, R:= sup \\{<P,dt(t>,V(l),uV^,v,Vv){t)\\^^,^y 
*e[o,T] 

Remark 4.4. By Assumption 14.21 all the functions we encounter in the 
previous statement are in C^{[0,T]; (O)) . Yet, one can verify that the 
estimate is valid whenever its two sides are well defined. 

Theorem 14.31 is not enough for the purpose of proving high frequency 
estimates independent of k . We will need the following version. 

Theorem 4.5. The statement of Theorem \4-^ remains valid if the sys- 
tem (|4.1() is replaced with 

' g,{(P){dtp + v-Vp) + ^divv-'^dw{/3i{cl))Ve)=fu 

(^•^^ ' g2i(^){dtv + v-Vv) + ^Vp- fif32i<l))Av- fif3li^)Vdivv = f2, 

. 9-M{dt0 + v-Ve)+ G{cl),Vcl)) -v + divv- kI3^{<I))M = /a, 
where G is smooth in its arguments with values in . 

The proof of Theorem 14.31 relies upon two estimates. Because they 
may be useful in other circumstances, we give separate statements in ^4.31 . 
The proof of Theorem 14.51 follows from a close inspection of the proof of 
Theorem 14.31 and so will be omitted. 

4.2. Example. Some important features of the proof of Theorem 14.31 can 
be revealed by analyzing the following simplified system: 



(4.5) 



dtp+ -divu - -M = 0, 

e e 

dtv ^ — Vp = 0, 

e 

dt9 + diYv-f3/\e = {). 
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For the sake of notational simplicity we abandon the tildes in this subsection. 
Parallel to (|4.2j) . we suppose 

(4.6) p>l. 

To symmetrize the large terms in , we introduce Ve := v — V6 . This 
change of variables transforms 1)4. 5|) into 



dtp H — div Ve 

£ 



0, 



dtVe + - Vp - V div Ve + {P- 1) VA6I = 0, 

dte + dwve- ip -i)Ae = 0. 

We take the scalar product of the first [resp. the second] equation with 
p [resp. Ve]- It yields: 

(4.7) ~ \\{p,Ve)\\l2 + \\divVe\\l2 -iP- l){Ae, divVe) = 0. 

We take the scalar product of the third equation with —rjAO , where r] 
is a positive constant to be determined later on. It yields 

~\\{p,Ve,V^yO)\\% 

+ lldivT;^!!^^ _ (/3 _ 1 + ^)( A0, diYVe ) + v{P- 1) IIA^Hia = 0. 
Set rj := f3 — 1 , which is positive thanks to assumption 1)4. 6|) . We get 



~\\{p,Ve,^/p^Ve)\\l, + \\d\y^. 



iP-l)A9\\l,=0. 



Let t ^ . Integrating the previous identity from to t , and using the 
triangle inequality to replace Ve by f , we get 

ft 



(4.8) \\ip,v,ve)mh+ \\divv-pAe\\i,dT^Kp\\{p,v,Vem\\i,. 

Jo 

Here and below, Kf^ denotes a generic constant which depends only on /? . 

We thus have proved an estimate independent of e . To go beyond 
and obtain smoothing effect on divv it is sufficient to estimate A9 inde- 
pendently. The strategy is to incorporate the troublesome term divv [in 
the equation for 6] into a skew-symmetric operator. 

To do so introduce 



C := £Pp — and Vg := £v. 



We compute 



dtC + 



div Vr = 0, 



dtVe + -^vc + ^ve = 0, 

P£ p£ 

dte + -dwv,-PA9 = 0. 

£ 
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We now use the essential feature of the system, that is assumption (|4.6|) . 
Multiply the first equation [resp. the second] by — 1) [resp. to 

put the penalization operator in symmetric form. The energy estimate thus 
reads 

~ ||(Vi/(/3-i)C, v^^e, e)\\l, + /? iiv^iii, = 0. 

Integrate this inequality, to obtain 

(4.9) \\{Q,ve,e){ml^+ f\m\i, dT^Kf,\\{c,ve,9)mi,. 

Jo 

Since the coefficients are constants, the same estimate holds true for the first 
order derivatives of C ; and 6 . Namely, one has 

(4.10) \\V{C,v„d){t)\\l,+ fwV^dWl, dT^Kp\\V{(:,Ve,9)m\l,. 

Jo 

On applying the triangle inequality, one can replaced C with ep in the 
previous estimates. Finally, by combining (|4.8|) . (|4.9j) and (|4.1()j) . we get 

\\{p,v,9){t)\\^, + \\V{e,ep,ev){t)\\^, + \\divv\\l, + \\V9\\l^ dr^ ' 
(4.11) 

^ Ki, \\{p,v,e)m\L^ + ||V(0,ep,et;)(O)||^. . 

It is worth remarking that, to establish ()4.11|) . it is enough to prove (|4.7() 
and (|4.1())) . Let us now compare (|4.1H) with the estimate given in Theo- 
rem I4.;^l We see that the previous study fails to convey one feature of the 
proof, namely the usefulness of the additional smoothing effect for Vp . The 
reason is that we worked with constant coefficients. In the general case, an 
estimate for ||Vp||^2 dr is needed to control the left hand side of (|4.10j) 
(see Lemma [4. 121 below) . 

Let us explain how to estimate Jq* || Vp||^2 dr . Multiply the second equa- 
tion in (|4.5|1 by eVp and integrate over the strip [0, t] x D , to obtain 
t i-t 

\\Vp\\l, dr = - {dtv,eVp)dT. 
Jo 

Integrating by parts both in space and time yields 

[\edtv ,Vp)dT = - [\v, edtVp )dT + e[{v{T), Vp(r) )] ^I* 
JO Jo 

t-t _ 

= / {divv,£dtp)dT + £[{v{T),Vp{T))]lzl- 
Jo 

In view of the first equation in 1)4. 5(1 , we are left with 

llVplli^ dr = f WdiYvWl, -{d\vv,M)dT -e[{v{T),Vp{T) 
Jo 

All the terms that appear in the previous identity have been estimated 
previously. As a consequence the estimate (|4.11|1 holds true if we include 
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/o 11^^*1112 dr in its left hand side. By so doing, we obtain the exact analogue 
of the estimate given in Theorem 14.31 (for /i = and k = 1 ) . 

4.3. estimates. Guided by the previous example, we want to prove 
estimates for (p, v) and (0, ep, ev) . The strategy for proving both estimates 
is the same: transform the system (|4.H) so as to obtain L? estimates uniform 
in e by a simple integration by parts in which the large terms in 1 /e cancel 
out. Namely, we want to obtain systems having the form 

(4.12) Ci{v, <^)U - CiilJi, (t>)U + -^S{cP)U = F, 

where >Ci(f , — £2(/U, k, 0) is a mixed hyperbolic/parabolic system of equa- 
tions: 

LQ{(j))dt+ ^ Lj{v,(j))dj - ^ Ljk{fj.,K,(l))djdk, 

and the singular perturbation S{(f)) is a differential operator in the space 
variable which is skew-symmetric (with not necessarily constant coefficients). 

We first prove an estimate parallel to (|4.9|) . 
Proposition 4.6. Using the same notations as in Theorem \4-'% we have 

sup \\{e,ep,ev){t)\\^2 + ( I K\\Ve\\\2 + fi\\eVv\\\2 dt] 

(GfO.Tl \Jq J 



(4.13) 



^ C(iio)e^^(«) \\{e,ep,evm\\L^+C{R)(^j^ B,{f,U)dt^ 



\ 1/2 



where 

(4.14) B,{f,U) := ||(/3,e/i,e/2)IL2 \\{d,ep,ev)\\^2 



1/2 



Jo 



Corollary 4.7. The estimate (|4.13jl holds true with (/g Bs{f,U)dt 
replaced by 

e||(/l,/2)IL2 + 11/311^2 dt. 



Proof of Corollary \4. 7| given Proposition |^.6] . One has, for all A > 0, 
[ B,if,U)dt] sup mep,ev)\\^2 + X [ ||(/3,e/i,e/2)IL2 dt. 

Jo / ^ te[0,T] Jo 

Using 1)4. 13() we obtain the desired result by taking A large enough. □ 

Notation 4.8. Within the proofs of Proposition 14.61 and 14.101 Rather 
than writing C{R) and C(i?o) in full, we will use the following abbrevi- 
ations: C = C{R) and Co = C{Rq) . So that C denotes generic constants 
which depend only on R := || (</>, V0, z/V^(/), f , Vw) ^j^P^ and Co 

denotes generic constants which depend only on Rq := \\(t>{0)\\i^oo(py As 
usual, the values of Co and C may vary from relation to relation. 
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Proof of Proposition \4 b\ . All the computations given below are meaningful 
since we concentrate on regular solutions of (|4.1|) . Introduce 

C := egiPsp- gsPiO, % := ev and U := (^,Ve,9)\ 

where we simply write gi , [5i instead of (7i(0) , I3i{4>) . We first show that U 
solves a system of the form (|4.12|) . Expressing p in terms of C, and 6 , 
replacing v by e~^Vc and performing a little algebra, yields 

dtC + v-VC + ^^^d\y% = f[, 

92{dtVe + V ■ VVe) + \^ {nZ) + ^V(72^) " ^Ji|32^Ve " ^l|3lV d\wv, = e/2, 
g^[dte + vVe) + -^ diwv, - kP^M = /3, 

, 1 53/31 , 

where 71 := — — , 72 := — — and 

f'l ■■= eMi - Pif3 + KpsVPi ■ ve 

+ ep{dt{gip3) + V ■ V{gip3)) - 0{dt{g3pi) + v ■ (V^s/^i)). 

In order to symmetrize these equations we make use of the structural as- 
sumption (|4.2j) . that is /Ss — /3i > 0. Multiply the first equation [resp. the 
third] by 71 (/Ss — [resp. 72]. By so doing, we obtain that U solves 

(4.15) Ci{v, (t))U - £2(/", (t>)U + -^S{(P)U = F, 

with F:= (7i(/33-/?i)-V{, e/2, 72/3)^, A(t',</') := Lo(m + E ^iK 
where 



71 







-^^0 := I g2ld I and Lj := vjIq. 

7253; 

Moreover, the singular and viscous perturbations are defined by 

/ 7idiv \ /O 

5= V(7i-) V(72-) ,-^2= ;U/32A + ^/3fvdiv 
\ 72div / \0 K72/33Ay 

The end of the proof proceeds by multiplying by U and integration by 
parts. The point is that the large terms in cancel out for S is skew- 
symmetric. Let t E (0, T] . If we multiply (|4.15|) by U and integrate from 
to t , we obtain 

(4.16) (Lo(0)Z?, U){t) - 2 [\c2{4>)li, U) dr = {Lq{<P)U ,U){Q) + Y{t), 

Jo 
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where ( • , • ) denotes the scalar product in L^(D) and 

Y{t):= n[dtLo{^)+ ^ djLj{v,(^)}u,ll)dT + 2 f\F,U)dt. 

Here we used the symmetry of the matrices Lq and Lj . 

We begin by estimating the right-hand side. One easily gathers 

dtLo{4>) + ^ djLj{v,(t)) 

(where C is as in Notation 14.8(1 and 

(4.17) \{F,U)\ ^ \\F\\L2\mL2^C\\{f[,ef2,h)\\L2 

Furthermore, since by definition /( is a linear combination of e/i , /s , ep, 
6 and kV^ whose coefficients are estimated in norm by a constant 
depending only on R, we get 

(4.18) \\f[\\^C\\{h,eh)\\^,+CmemL^+C\\KVe\y. 
Moreover, observe that 

(4.19) 0,£p,£vf = M{(j))U, 

where M{(j)) is an nx n invertible matrix. As a consequence, there exists a 
constant C such that C'^ < II (6', ep,^) 11^2 ^ C||Z^||^2. Hence, the 

estimates (|4.17() and ()4.18|) result in 

(4.20) K z?)K CB,{f, u) + c ml, + c WkVOW^, im^, , 

where B^{f, U) is given by (j4.14|) . For all A ^ 1 , the last term in the right 
hand-side of the previous estimate is bounded by 

AC||^?||i. + y llV^lli.. 

Consequently the right-hand side of (j4.16|) is less than 
(4.21) 

C7o||Z?(0)||i2+AC /" \\U\\l,dT + ^f \\Ve\\l,dT + c[ B,{f,U)dr. 
Jo ^ Jo Jo 

We want to estimate now the left-hand side of (|4.16|) . In this regard, the 
most direct estimates show that 

{Lo{<P)U,U)^\\Lo{(t>)-'\\li 

Similarly, using Lemma l4 . 91 b elow . one has 

-{C2{(t))U, U) ^ Kiim llVv^f + Knm ||V^f - C \\U\\\2 , 

where m := min{ H/Jg^^ ||^i , ||(/32 + PIY'^W'lL , WPs'^Wl^} and K is a generic 
constant depending only on the dimension. 
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To estimate the L°° norm of Lq ""^ , ^ ~^ (^2) ^ ^'^^ ^ > write 



(4.22) ^Co + t sup ||F'(0(r))||^^ l|5t</.(r)||i^ 

re[o,t] 

Consequently, the left hand-side of (|4.16p is greater than 

(4.23) Coe-^^{||Z?(t)||i2 +^Vl|Viyef + K||V^f dr} -C^*||Z?||i2 dr. 

From this and the assumption k ^ 1 , we see that one can absorb the third 
term in ()4.2ip by taking A large enough. Then, Gronwall's Lemma implies 

\\^^{t)\\h+ [\\\Vvef + K\\Vef dT^Coe^^\\Um\l, + C f B,{f,U)dT. 
Jo Jo 

We next claim that the previous estimate holds true with U replaced 
by {6,ep,ev) . Indeed, using the identity (|4.19|) and the bound ()4.22() . one 
has ||(6',epi, e'u)||^2 ^ Coe^^ W^Wl^ ■ To complete the proof, take the square 
root of the inequality thus obtained and take the supremum over t G [0, T] . 

□ 

Let us briefly recall a standard estimate used in the previous proof. 

Lemma 4.9. There exists two constants Ki = Ki{d) such that for all Lip- 
schitz functions C = C{x) and rj = r]{x) such that C > and + C > 0; 
and for all vector field u £ H^{0) , there holds 

2 K2M'^ ,,2 

(4.24) — (CAti + r^Vdivu, u)h-ixH^ ^ Kim \\Vu\\j^2 IpIIl2 > 

where m = mfx(zo{C{x), ((x) + ri{x)} and M = ||VC||ioo + ||Vt/||j;^oo • 

Proof. Decompose u as u + Vcp := Qu + (/ — Q)u where Q is the Leray 
projector onto divergence free vector field ^, so that 

(An + r?V div n = CAu + (C + r?) VA0. 

It results from divn = that 

-(CAn + r/Vdivti, u)H-ixm = (CVu, Vu) + { {C + rj) Ac/) , A^) + R, 

where R is such that \R\ ^ K^M \\'^u\\^2 \\u\\j^2 ■ 

To infer the desired bound, we use two usual inequalities. For all A ^ 1 , 
one has \R\ ^ (m/X) || Vu||^2 + {XK^AP /m) ||u||^2 • We next recaU the most 
simplest of all C alder 6n-Zygmund estimates: ||V^(/>||2,2 ^ K4 \\A(j)\\j^2 ■ □ 

We next prove an estimate parallel to (|4.7|) . 



^This trick is due to R. Danchin. 

24 



Proposition 4.10. Using the same notations as in Theorem we have 
sup \\{p,v){t)\\^2 + [ / fi\\Vv\\^2 + K\\divv\\j^2 dt] 



^ C(i?o)e^^(^) ||(p,?^)(0)||^. +C(i?) f \\{h.h)\\L-^^\\h\\n^ dt 

1/2 







+ C(i?o)e^^(^)| sup + ( /^k||V0||^i dt 

Ue[o,T] " V^o 

where v := ^/JT+~k and = \\-\\x^2 + v W^i ■ 

Remark 4.11. From now on, we make intensive and implicit uses of the 
assumptions /i ^ 1 and k ^ 1 . In particular we freely use obvious estimates 
like /i ^ yjjl ^ V . Similarly, we use the estimate v \\^u\\i^2 ^ ll^llj^i without 
mentioning it explicitly in the proofs. 

Proof. Introduce Ve := v — nPiVO. Performing a little algebra we find 
that V := {p,VeY satisfies 

;^j(,v.„.vv).l(« t)v-(M)-^, 

where 

£22 := Kg2piV{g^^ div • ) + /^/32A + /u/3^Vdiv, 
and F= (/i,/2 + /^)* with 
(4.25) 

:= ^lKP2^{Pl^0) + fiKP^V dWiPiV 9) - Kg2{dtl3i + v ■ Vpi)Ve 
+ «;52/5iV(k53~' dwiPiVd) - Kg^'PsAd- g^' h) + Kff2/3i(Vi;)*V^. 

The thing of greatest interest here is that —£22 is a differential operator 
whose leading symbol is greater than Ck^* • + Cji |^|^ . 

As before, the proof proceeds by multiplying by V and integrating on the 
strip [0, t] X D . Then the analysis establishing that the right [resp. left] hand 
side of (|4.16j) is smaller than (|4.21|) [resp. greater than (|4.23p ] also gives: 



(4.26) 



II'>^(*)IIl2 + / M l|V'Ue|li2 + K ||divt!e||^2 dr 
JO 



^Coe^^||V(0)||i2 + C f\\V\\l2 dr + C f\{F,V)\dT. 

Jo Jo 



Let us estimate the last term in the right-hand side of ()4.26|) . Using 1)4. 25(1 
one can decompose /2 as /2 1 + \/^^/2 2 where 

||/2,ilL2 ^ Ck wMhi + Ck \\ve\\L2 + cv^ifi + k) \\v^\l2 , 
||/2,2lL2 ^ Ck \\ve\\^2 + cMfi + 1,) \\v^\^2 . 
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By notation, ( F , V ) = ( /i , p ) + ( /2 , ) + ( /£i , ) - ( /g^s ' \/^div ) . 
We thus get, for all A ^ 1 , 

\{F,V)\dT^XC{ I ||(/i,/2)IL2+k||/3||^^i dr' 







If 1 ~ 

(4.27) +-/ K\\dw Ve\\l2 dr + - sup ||V(t)||^2 

A Jo A T-e[o,t] 

+ XC [ \\{Kve,^/^{^l + K)v^e\\l2 dr. 

Jo 

Replacing with v in H4.26() - ()4.27|) . taking A large enough and applying 
Gronwall's lemma leads to the expected result. □ 

4.4. End of the proof of Theorem 14. 3L From now on, we consider a 
time < T ^ 1 , a fixed triple of parameter a = (e, £ A and a solution 
U = {p,v, 9) of the system ()4.1() . We denote by Rq and R the norms defined 
in the statement of Theorem 14.31 (see (|4.3() ). We also set v := + k . 
Introduce the functions N,w,y,Y,Z: [0, T] — > M_|_ given by 

(4.28) Ar(t):=||^||^„^^^, 
and 

/ft ^2 ~2 V^^ 

w{t) := sup \\{p,v){t)\\^2 + [ / ^llV^ll^a +K||divv||^2 dr , 
re[o,t] \Jo J 

~ ^ ^ / /■* ~ ^ 

y(t):= sup ||(^,ep,et;)(r)||^2 + / K\Vef^2 ^ ii\\eVvf^2 dr] , 

y(t) := sup ||V(^,ep^,e£^)(T)||^2 + / k||V20||2, +^||eV2?y||2, dr , 

refO,tl \J0 / 



1/2 

^(i):=( / llVplL^dr" 



Note that A^(T) ^ ?x;(T)+z^z(T)+7/(T)+z^y(T) . Since we have estimated 
wi^T^ and y(T) , it remains only to estimate vYiT^ and i'z{T) . Parallel 
to ()4.1U|1 , we begin by establishing an estimate for Y . 

Lemma 4.12. There exists a constant Cq depending only on Rq and a 
constant C depending only on R such that for all t G [0, T] and all X ^ 1 , 



(4.29) 



Y{tf ^ Coe^^Y{Of + XC f Yirf dr + ^ [ 

Jo Jo 

+ [c l^\\h\\Hi + \\(.fij2)\\Hi dr^ . 



\\{divv,Vp)\\l2 dr 
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Proof. The proof, if tedious, is elementary. Indeed, the strategy for prov- 
ing (|4.29j) consists in differentiating the system (|4.1|) so as to apply Propo- 
sition with {6,ep,ev) replaced by V{6,ep,ev) . 
Let 1 ^ j ^ d and set 

Uj := {pj,Vj,ej) := {djp,djv,dj6). 

We commute^ the ith equation in (|4.1|) with gi{djg~^-) . It yields 

9i{^){dtPj + V ■ Vpj) + ^ diYVj - div(/3i((/))V%) = Fi, 

g2{(t>){dtVj + V ■ Vvj) + ^Vpj - pP2{^)^Vj - ///3^((/>)V divvj = F2, 
^ g^{(t>){dtej + V ■ VOj) + dwvj - Kp3(.^)Aej = F3, 
where, for i = 1,2,3, the source term Fi is given by 

Fi := Fl + Ff := gid,{gr^ h) + giPf , 

with 

Fl := - djv ■ Vp- ^djg^^ divv+ '^djg^^ div(/3iV^) + "^g^^ div{dj(3iV9), 

F| := - djV ■ Vv- ^djg^^\/p + pdj {g^ ^ P2)^v + pdjig^^ l3l)V diyv, 
F| := - djV ■ V6- djg^^ divv + Kdj{g^^ 13^)^. 
Proposition 14.61 implies that 

(4.30) Y(tf ^ Coe^^ \\Um\l^ + ^ f ^^(^' ^i) 

where U := {V6,£Vp,eVv) and is defined by ()4.14p . We now have to 
estimate {F, Uj ) . The source terms are directly estimated by 

\\{Fi,eFleFi)\\^,^C\\h\\H.+C\\{h,f2)\\Hi, 

\\{Fi,eFleFi)\\^,^C\\U\\L, + C\\{dWv,VmL^+C\\{epV^v,KV^9^^^ 
The last estimate implies that, for all positive A and A' , one has 
B,{F\U,) ^ (l + A + A')C||^^||2, 

+ ^ IKdiviJ, VfOlli, + ^ \\{epV'v,t,V%\\l, . 

Moreover, the term Jq\Bi;{F^ ^ C/j)| dr is estimated as in the proof of Corol- 
lary U7J With these estimates in hands, we find that the last term in the 



■^Which means that we commute dj with the ith equation (i = 1,2,3) premultiphed 
by ^iid we next multiply the result by Qi . 
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right-side of (|4.3Uj) is less then 

[\i+x+x')c\\u\\i, dT+^ sup ml, 

Jo ^ r6[0,t] 

(4.31) + l^^\\(divv,Vp}\\l^ dT+^ l\\{ef,V'v,KV^\\l, dr 

+ >^"(c j\\h\\m + \\{h,f2)\\Hi dr^. 

Yet, by definition, sup,g[o,t] ll^lll^ + /o IK^^/^V^iJ, kV^^)!!^^ dr ^ Y{t)\ 
Hence, taking A' = A" large enough, one can absorb the second and fourth 
terms in (|4.31() in the left hand side of H4.30() . thereby obtaining the desired 
estimate □ 

To obtain a closed set of inequalities it remains to estimate z{t) . 

Lemma 4.13. There exists a constant Cq depending only on Rq and a 
constant C depending only on R such that for all t G [0, T] , 

(4.32) u^z{tf ^ Coe''''{u'Y{tf + w{tf} + C(^I^ \\{fij2)\\L2 dt^ . 

Proof. Set := [0,t]x]I]) and denote by (( , )) the scalar product in L^(Qt) . 
Multiplying the second equation in (|4.1I) by ei^^(72^^Vp, and integrating 
over Qf yields 

j^^iig^'Vp, Vp)) = 

(4.33) - ei^^ii (dt + v V)v, Vp)) + e^i^^Hg^^ (32Av , Vp)) 
+ efiu^{{g^'plVdivv, Vp)) + eu^{{g^' f2,Vp)). 

The most direct estimates show that 

5^^/52 Aw, Vp)) ^ ||5^^/32|li.c(Q^) \\efiuV'^v\\^2^^^) IW^Ph^nt) ' 
efii^^{{g2^4vdiyv, Vp)) ^ {{92'^ (^Ih-^iUt) l^^^^^ L^Qt) IW'^Ph^in,) , 
£^'^{{92^f2,Vp)) ^ \\92^\\L--{nt) II/2IILIL2 lk^^Vp||^o^2 , 
where we used the shorthand notations 

II • ll^ox := sup II • 11^ and || • ||^p^ := ( [ \\ ■ ||^ dr^ . 

* Te[o,i] * \Jo J 

With the bound (|4.22|) . the previous four estimates and (|4.33p imply 
(4.34) 

ll^Vplli^Ca) ^ ^oe^^ Vv^k idt^v- V)v, Vp))\ 

+ Coe^^(||e^i/V2?;||^2(f,^) + II/2IL1L2 + \\eu^Vp\\cOL2y . 
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Note that the second term in the right-hand side of the previous estimate is 
bounded by the right-hand side of (|4..S2j) . It thus remains to estimate the 
first term. To do so integrate by parts in both the space and time variables 
to obtain 



ev 



2/ 



\{dt + vV)v,vp)) 



= eu^ ((div v,{dt + v- V)p)) 

-eiy'^{dwv{t),p{t))+eiy'^{divv{0),p{0)) 
+ ev^ {{ Vv , Vp (8) w )) — eu'^ {{ div v,v- Vp )) . 

The last two terms in the right-hand side of H4.35() are estimated by 

Kt ||Vv||^Oioo \\vfcOL2 + \\ei^Vp\\ljOL2 < TCw{tf + v'^Y{tf, 

and the second and third terms are estimated by 

WevVvWl,^, + WpWI,^, ^ u^Y{tf + w{tf. 

In particular, the sum of the last four terms in the right-hand side of (|4.35|) 
is estimated by e^^ {w{t)'^ + z/^y(t)^} . This brings us to estimate the first 
term. Using the first equation in (|4.1() . we get 

ei^^((div?;, {dt + v ■ V)p)) = — i^^((divS', g^^ divw)) 
+ u^K{{d\Yv, <7f ^ div(/3iV^) )) + eu^idiYv, g^^ h ))■ 
Then, the analysis establishing 1)4. 34() also gives 
ez^2|((divv, {dt + v-V)p))\ 

^ C,e^^[v^ ||(div?^,KV2^)||i.(^^) + v^K\\{Ve,eVv)\\loL, + H/iHii^^}, 

which in turn is bounded by the right-hand side of l|4.32() . The proof is 
complete. □ 

We thus have proved a closed system of estimates. Indeed, by taking 
appropriate combinations of the previous estimates, we see that there exists 
a constant Cq depending only on Rq and a constant C depending only on 
R such that for all t G [0, T] and all A ^ 1 , the norm N{t) (as defined 
in H4.28|l l satisfies the estimate 

N{tf ^ Coe^^N{<df + XC [ iV(r)2 dr + ^N{tf + CF{Tf, 

Jo ^ 

where F{T) := J^* ||(/i, /a)!!^,!^ + WMhi dt . 

Again, take A large enough and apply the Gronwall's lemma. The proof 
of Theorem 14.31 is complete since ||f^||-^o(2") ^ N{T) . 
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5. High frequency regime 



To establish Theorem 12.71 the crucial part consists in obtaining a priori 
estimates in Sobolev norms independent of a G A . Theorem l4.3l Drovides the 
basic estimates. However, the estimates of the derivatives also require a 
careful analysis. Indeed, the classical approach, which consists in differenti- 
ating the equations, certainly fails since it reveals unbounded terms in . 
However, one can follow this strategy in the high frequency regime where 
the parabolic behavior prevails. 

More precisely, given a smooth solution U = {p,v,9) to (|2.1j) . we will 
estimate the Ti.^{T) norm of (/ — Jh)U , where TC'^{T) is as defined in Def- 
inition 12.41 and, in keeping with the notations of ^ {Jh\h£ [0,1]} is a 
Friedrichs' mollifier and A'* := (/ — AY^"^ . In order to make our energy 
estimates applicable, the main difficulty is to verify that the commutator of 
(/— J/j)A® and the equations <\2.1^ can be seen as a source term. To do so we 
first note that for h = 0{e) , one can gain an extra factor e in the commuta- 
tor estimates [see (|5.2|) below]. Yet, this costs a derivative. To compensate 
this loss of derivative, we use in an essential way the parabolic behavior of 
the equations. Consequently, we search h under the form c{fj,, k)£ . Since 
for our purposes the main smoothing effect concerns the penalized terms 
divi; and Vp, we take c{^,k) = ^ k. 

Proposition 5.1. Given an integer s > 1 + d/2, there exists a continuous 
non- decreasing function C such that for all a = {e, fj,,K) € A, all T G [0, 1] 
and all U = {p,v,9) G C\[0,T]; H^{B)) satisfying (tTTl) . 

(5.1) 11(1 - J..)C/||^.(r) ^ C(J7o)e^^(^), 

where u := ^/jr+~K , Qq := \\U {0)\\')^b ^ , Q := 11^^117.^8 (j^) (see Definition \^^ ) . 

Remark 5.2. Note that we establish estimates for the exact solutions 
of (|2.1|) and we do not estimate the solutions of the linearized system 1)4. 1|) . 

5.1. Preliminaries. To avoid interruptions of the proofs later on, we now 
collect a few nonlinear estimates we use throughout this section. 
Recall that ||-||j:^<t := ||-||j|^<T-i + Q \\-\\fjcr ■ 

Lemma 5.3. Let s > 1 + d/2 . There exists a constant K such that for all 
/i G [0, 1] , a// £» ^ , all f G iJ''+^(D) and all u G H^B) , one has 

(5.2) \\[f,{I-Jh)A']u\\^, ^ {e+h)K{\\f\\^,,^\\u\\HB + \\f\\^B+,\\u\\^^}, 
where \\f\\^,,^ := + ||V/||^^ . 

Remark 5.4. This inequality provides a way to gain an extra factor h 
[with g = h] oi a derivative [with g = 1]. In this respect it is like the 
estimate (|3.11|) . to which we will return in a moment (see Lemma 16.6)1 . 

Proof. To prove this result we need a tame estimate version of ()3.1U() . We 
use the following result, which is Proposition 3. 6. A in [301 (see also [25]). 
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Let m > . For all Fourier multiplier Q with symbol q G S"* and for all 
cr ^ 0, we have 

(5.3) II [/, Q]it||j|^<T ^ 11/11 + 11/11 

where the constant K depends only on m , s , cr , d and a finite number of 
semi-norms of q in S"^ . 

Since the support of the symbol 1— jh is limited by |^| ^ the most 
direct estimates show that {h~^{l — j{h^)){S,)^ \h G [Oil]} is uniformly 
bounded in the symbol class S^~^^ . Thus, the commutator estimate H5.3|) . 
applied with (m, a) := (s + 1,0), implies 

(5.4) II [/, (I - Jh)A']u\\^, < h{\\f\\^^.^ \\u\\j,s + 11/11 J^s+l ||'U|| 2,00) . 

On the other hand, the family {(1 — j(/iO)(0* 1^ ^ [0)1]} is uniformly 
bounded in , so that 1)5. 3(1 applied with (m, a) = (s, 1) implies 

(5.5) II [/, (/ - JhW]u\\^, < 11/11^,,^ ||n||^. + 11/11 fjs+l \\u\\ j^oo ■ 
Combining ()5.4() with ()5.5() multiplied by g, yields 1)5.2(1 . □ 

We next prove two Moser-type estimates for the norms |M|t7<T . 

e 

Lemma 5.5. Let s > 1 + d/2 . There exists a constant K such that for all 
g ^ and for all Ui G Hg[ 



(5.6) ||nin2||frs ^ K \\ui\\rTs \\'U2\\u!> ■ 

This result extends to vector valued functions. 

Proof. Using the standard tame estimate for products at orders s — 1 and s , 
we get 

||^^1^^2|liys-i < ll^illLoo ||^i2|li^s-i + llttlllns-i ||^i2|lL°o ) 
ll'"l'"2|liys < ll^^llli^oo \\u2\\hs + 11^^211^00 . 

Using the definition ||-||j:^s := ||-||j|^s-i ||-||j|^s , and putting the parameter g 
in appropriate spots yields 



Since s— 1 > d/2, the Sobolev Theorem implies ||n||^oo < < H^^H/fs • 

Which completes the proof. □ 



Lemma 5.6. Let s > 1 + d/2 and F : M" —^Cbea C°° function such 
that F{0) = . Then, for all g ^ and all u £ L[^{D) with values in M" , 

(5.7) \\F{u)\\j,.^^C{\\u\\hs), 

where C(-) depends only on a finite number of semi-norms of F in C°° . 
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Proof. Recall that, for all ctq > d/2 and all C°° function F satisfying 
F(0) = 0, 

(5.8) \\F{u)\\h.,^C{\\u\\^^)\\u\\^^,. 
Using this estimate at orders s — 1 and s , one has 

(5.9) \\F{u)\\^s^C{\\u\\^^)\\u\\hs. 

This in turn implies the desired estimate. □ 

Lemma 5.7. Let s > 1 + d/2, F: ^ C be a C°° function such 
that -F(O) = and Q he a Fourier multiplier with symbol q £ S'^ . For 
all vector-valued function u £ H^{0) one has 

(5.10) \\Q{F{u)) - F'{u)Qu\\^, ^ C{\\u\\j,.), 

where F' is the differential of F and C(-) is a smooth non- decreasing func- 
tion depending only on s , d, a finite number of semi-norms of q in and 
a finite number of semi-norms of F in . 

Proof. To establish (|5.10() we use the para-differential calculus of Bony j^j. 
Denote by Tf the operator of para-multiplication by / . Starting from 

F{u) = Tp,^^)U + R{u; x, Dx)u, 

where R{u; x, Dx) is a smoothing operator (see (|5.1H) V we obtain 

Q(F(n)) -F'{u)Qu = {Tf,^^^-F'{u))Qu+ [Q,Tf,^^)]u + QR{u; x, Dx)u. 

The claim then follows from the bounds (|3.16j) - (|3.18|) and the estimate 

(5.11) \\Riu;x,Dx)u\\H,+^ ^ C(||n||^,). 

See dZl Th. 10.3.1] and [13 Corollary 9.3.6] for the proof of □ 

5.2. Localization in the high frequency region. To simplify the pre- 
sentation, we fix a real number s strictly greater than 1 -|- d/2 . 
To proceed further, we need some more terminology. 

Notation 5.8. For ah h G [0, 2] , define 

Qh := {I - Jh)A'. 

Hereafter, the parameter /i is a product eu with (e, ly) G [0, 1] x [0, 2] . 

Litroduce next the commutator of the equations 1)2. 1|) and Qsu ■ 
Notation 5.9. Given a = {e, fj,, k) £ A, u e [0,2] and U = {p,v,9), set 

/r:HF(^) ■■= [9i{<t>),Qe.]dtP+[gi{(t>)v,Qsu]-Vp -'^[Bi{4>),Qsu]9, 

f2:RF(U) ■■= [52(0), Qe.]dtV + [g2{(p)v, Q,,] -Vv -fl [B2{^), Qeu]v, 

f^^YiU) ■■= h(0), Qeu]dte + [g^v ■ Vd; Q,^ - k [B^i^), Qeu]0, 
where i;^ is a shorthand notation for {9,£p) and 

{g3{(l))v ■ Ve- Qeu] := g3{(l})v ■ VQeud + 53(0)Ve • QeuV - Qeu{g3{<t')v ■ VB). 
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Remark 5.10. For the purpose of proving estimates independent of k, 
we do not consider the exact commutator of the third equation in (|2.1j) 
with Q^^ . 

Let us recall that 

5i(0) :=XiMdiv(/?(e)V-), 

B2{(t>) ■■= X2{ep) div(2C(e)I?-) + X2{ep)V{m div •), 

Sm :=X3Mdiv(/3(e)V-). 

The following lemma shows that /^p can be seen as a source term. 

Lemma 5.11. There exists a smooth non- decreasing function C = C(-) 
such that for all a £ A , all T ^ , all v £ [0, 2] and all vector valued 
function U := {p,v,9) £ C\[0,T]; H'=^ (B)) , 

||/r;HF(f^)Li ^C{R){l + \\edtp\\Hs+K\\e\\Hs+2}, 

\\f2:uFiU)\\Hi ^C{R){l + \\edtv\\H + f,\\ev\\Hs+.}, 

\\ftnAu)\\^, ^c{R){i + \\dte\\Hs +K 11^11^.+,}, 

where R := + H^H 

Remark 5.12. We will apply this lemma with i/ := ^/i + k . Yet, we prove 
estimates independent of G [0, 2] to explain why the frequency space is 
cut around |C| ~ l/{e^J + k) . 

Proof. We make intensive use of the following obvious observations. Firstly, 
using the Sobolev embedding Theorem and the very definition of the norms 
IMI o-a , one has 

(5.12) ||u||^oo < ^ Ikll/fg and ||n||^yi,oo < \\u\\fjs < . 
Since e ^ 1 and 1/^2, directly from the definition of IMIh-ct , one has 

(5.13) u \\u\\j^a ^ ll^^llji/CT ^ 3 ll^llji^CT and ||en||^s+i ^ ||m||^s+i . 

Note the following corollary of the second inequality: ||(ep, e7;)||j:^s+i ^ R. 
This in turn implies := ||(0,ep)||j|^s+i ^ R. 

Step 1: Estimate for /"hf(U)- a) We begin by proving that 

(5.14) II [5i(</>), Qeu]dtp\\Hi^ ^ C{R) WedtpWfj. • 
Applying the commutator estimate ()5.2() with h = g = ei^ , we have 

\\[9i{(l^):Qeu]dtp\\jji^ <ejy\\gi{4>)\\w^,o. WdtpWns + £^\\9i{(I>)\\h^+i WdtpWio^ 

^ Il5l (</>)! I \\^dtp\\Hs + ||5l((A)ll/f=+i W^dtpWLoo , 

where gi is defined by gi = gi — gi{0) . The estimates (|5.12|) imply that 
the right side is less than ||^i(i;^)||^s+i ||e9(p||j|^s • Using Lemma f5.61 we 
obtain ((/>)! ^ C j^s+i) ^ c\r) . This proves (f5Jl|) . 
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b) Next, we prove that 

(5.15) II [gi{^)v, Qeu] ■ ^vWhi^ ^ 

Again, this follows from the commutator estimate ()5.2() applied with h = 
Q = ev . Indeed, it yields 

II [9i(0)^^, Qei^] -VplLi < \\h{4>)v\\w^,oo\\eVp\\^sM\£gi{(l))v\\jjs+i\\Vp\\^oc, . 

Ill- -iiiJJgjy }j u 

The first term in the right-hand side is estimated by C{K) since 

Il5i(</')^^llvyi.°° ^ C'(ll<Alliyi,°= > Iklli4^i.'-) ^ C'dl'/'llns , ^ C{R), 

and ||eVp||j|^s ^ 11^11//'=+^ ^ 

It thus remains to estimate ||e^i((/))z;||j:^s+i . To do so use Lemma 15.51 
and 15.61 to obtain 

lk5i(0)v|lH=+i ^ II5i('A)IIhj+i Ik^ll^j+i ^ c{R)R. 

c) Let us prove that 

(5.16) -\\[B,{(l)),Q,,]e\\^, ^C{R)\\Ke\\Hs+2. 

£.111- J ii-Uei/ 

By definition Bi{(f)) = Xi{^P) div(/3(0)V • ) , so that one can decompose the 
commutator \^Bi{(f)) , Qsy\6 as 

(5.17) [£i{(fy, Qeu\^e + [£2{(t^,V(fy, Qeu] ■ Ve, 

where fi((^) = xiiep)m and £2{^,V^) = xi(ep)V/3(0) = xi{ep)/3' {9)^9 . 
The commutator estimate ()5.2p [applied with h = g = ev] implies 

II [£i, Qeiy]^9\\^-^ < £ \\£i\\yyi,oa \\A9\\jjs + e ||<Si||j|^s+i ||A6'||^oo , 

1^1,00 II V0|| ji^s + e \\£2\\fjs+i II V0|| j^oo , 

where £i := £i — <?i(0) . We gather 

l|fi(<A)ll^i.oo ^C(||0||^.)^C^(^), IIA^II 

Hf, ~ II^IIhs+2 ) 

\\£i{4>)\\h^+^ ^ C{\\4>\\h^+.) ^ C{R), \\Ae\\^^ < \\0\\j,s+2 , 



\u£2{<p,v<p)\\w^,^ ^ cm\j,s+,) ^ c{R), \\ve\ 



H= ^ 11^11/^5+2 



\\£2{4>,^ct>)\\Hi+^ ^ C{R) \\e\\Hs+2 , liv^ii^.. < ii^ii^. ^ R. 

The first six inequalities follows from Lemma 15.61 the Sobolev Theorem 
and the estimates ()5.12|) - H5.13p . To estimate ||<?2(i?i*, Vi;^)||j|^s+i , we first use 
Lemma 15.51 to obtain 

(5.18) ||f2(<A, V<A)||^.+i < (1 + \\xi{ep)\\Ht+^){l + \Wm\m+^) l|V0||^2+i , 

and next use Lemma 15.61 to bound the i?^"*"^ -norms of xi{^9) and 13' {6) by 
C{R) . This yields the desired bound since 11X76*11^:^8+1 < ||^||j:^s+2 . 
The claim ()5.16p then easily follows from these estimates. 
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Step 2: Estimate for /s^hfI^) • Note that one can obtain /2hf(^) from 
fi'^p(U) by replacing p by w and 6 by ev . Therefore, we are back in the 
situation of the previous step, and hence conclude that 

\\[g2{^),Qsu]dtv\\^, ^C{R)\\edtv\\Hs, 

\\[92{(p)v,Qeu] -Vz^ll^i ^C{R), 

" " -'--'-£1/ 

II [B2{(p), Qeu]v\\rri ^ C{R) \\fiev\\fjs+2 . 

Step 3: Estimate for hf(^) • estimates for the first and the last 

terms in f^'^p{U) can be deduced by following the previous analysis. We 
only indicate the point at which the argument differs: we use the commu- 
tator estimate ()5.5() with h = eu and g = u (instead of g = ev). 
Let us concentrate on the second term. We claim that 

\\{g^{<t>)vVe-Q,,]\\^,iiC{R). 

We first decompose {g'i{4')v -yO^Qeu] as 

(5.19) [g:i{4>),Qeu\v-Vd + g^{<i))Z 

where Z := v ■ VQeyO + V9 ■ Q^uV - Qeu{v ■ VO) . 

Applying the commutator estimate 1)5. 2|) with h = eu and g = u , and 
using the estimates 1)5. 12() , we get 

II [gsi^), Qeu]v- ve\\^, < ||53(</')IIhj+i 11^ • V0||^. . 

Using Lemma 15.51 and 15.61 we find that the right-hand side of the previous 
estimate is dominated by C(||0||j:^s+i) ||w||j:^s ||V0||j:^s ^ C{R). 

Since ||5(3((/))Z||^i ^ 1153(0)11^/1.°° II^IIh^ ^ C'(i?) ||Z||j:^i , to control the 
Hl-norm of the second term in (|5.19)) . only the estimate of is miss- 

ing. We split ll^'ll^i as ||Z||^2 + , and we prove that ||^||^2 and 

are both estimated by C{R) . To estimate ||Z||^2 , note that 

(5.20) Z = Ve- QevV + {v, Q,ev\ ■ V0. 

The second term in the right-hand side is estimated by way of the commu- 
tator estimate p.lOp . Indeed, applying Lemma 13.41 with ctq = = s and 
(T = s — 1 , we get 

\\\v.Qev\-^Q\\L^ < \\v\\hs \\Ve\\Hs-i ^ \\v\\hs IMhs ^ R^. 

As regards the first term in the right-hand side of 1)5.20(1 . write 

\\Ve-QeuV\\L2 ^ IIV^II^oo \\QeuV\\L2 < \\0\\hs \\v\\jjs ^ R^ . 

Moving to the estimate of ||Z||^i , remark that 
uZ = F'{u)QeuU-Qe^{F{u)), 
with u = (ni,M2) := {vjvVO) and F{u) = uiU2- Hence, ()5.10() yields 

u\\z\\^, ^c{\\iv,i^ve)\\Hs) ^ciM^. + Wew^s^.) ^ciR). 
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The previous estimates imply that ||^||fri is controlled by C(R) . This 
completes the proof of Lemma 15.111 □ 

Remark 5.13. Let us explain the reason why we assume that /3 depends 
only on 6 . Had we worked instead with general coefficient (3 depend- 
ing also on ep, the corresponding inequality H5.18() would have involved 
||eVp||j|^s+i . The problem presents itself: lkp||j;^i(o T_ff°+^) is not controlled 
by the norm f , 9)\\-^^s (^x) ■ 1^ is possible to get around the previous prob- 
lem, yet we do not address this question. 

In view of Lemma [5. Ill we are led to estimate dttp where ip := {6, ep, ev) . 

Lemma 5.14. There exists a continuous non- decreasing function C(-) such 
that for all a £ A, all T ^ , all u £ [0,2] and all {p,v,0) G W^(T) 
solving 1)2. 1(1 . the function -0 := {9,ep,ev) satisfies 

(5.21) \\dti'\\Hs^C{R){l + R'}, 

where 

R:=\\{p,v)\\jjs+i + \\e\\jjs+i , 

R! := V ||(divf , Vp)||^s -|- ^/Ji\\Vv\\^s+l + a/k || V6'||^s+i . 

Remark 5.15. This estimate plays a key role for the purpose of proving 
estimates independent of /i and k. Indeed, this is the only step in which 
we use the additional smoothing effect for d\vv and Vp. More precisely, 
the fact that the estimate ()5.21|) is tame [linear in R'] allows us to control 

\\^M\l-^{0,T;HI) by C{\\ip,V^G)\\ni(T)) ^OI Sl\ V ^ ^fi + K. 

Proof. In order to establish 1)5. 21(1 . observe that 

(5.23) G(0) {dtiJ + vVij) + LU + B^^^{ct>)^ = 0, 

where ip and U are identified with (ep,ev,6Y and {p,v,9Y , respectively, 
and 

fgi 0\ /O div 0\ /O kBi 

G := [0 g2ld , L := V , B^,, := \ fxB^ 
\0 gsj \0 div 0/ \0 kB^^ 

Since max{/i, k} ^ 1 , one can easily verify that there exists a family 
{Fa I a £ A} uniformly bounded in C°° , such that Fa{0) = and 

with E := (-y, ip, VV', div v, Vp, e/iV^f , . 
Using the Moser-type estimate ()5.9() . we find 

PiVII^,. ^c(||H||^._0{i + H|s||^.}. 

To complete the proof note that, directly from the definitions (|5.22p and 
the assumption {/J,, k) G [0,1]^, one has ^ R and ^ 

R + R'. □ 

We are now prepared to prove Proposition 15.11 
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5.3. Proof of Proposition IFTTl Set v := ^/x + k ^ \pi and {p,v,9) := 
(Qei^P, Qeiy?^, Qey^*) • We first show that U := {p,v,6) satisfies (|4.4|) for 
suitable source terms /i , /2 , /a and coefficients /3i , /32 , ) /^s and G . 
It readily follows from Notation EH that {p,v,9) solves 

' g2{ct>){dtv + v-Vv) + ^Vp-fiB2{(P)v = f^'^^p{U), 

where = {d,ep). Set /3i := /32 := X2C, /?2 ^2??, /^s := Xa/? and 

G{(j),V(j)) = g3{(j))V6, to obtain that {p,v,9) satisfies ()4.4() where 

/i - fi:KFiU) + /{,HF + Qe^T^i With /{_HF := -^Vxi(ep) • {(3m0), 

h ■■= f2:HFiU) + /2,HF With /^,HF := ^X2M{2I?WC(e) + divWv{9)}, 

h ■■= ft^AU) + /3,HF + eQe.Ta with /^^hf := '^X3(ep)V/3(e) • V^, 

where Tj are as in system (|2.ip . 

By definition of the norm || -||^s (see Definition 12.4(1 . the Sobolev em- 
bedding Theorem and Lemma l5 . 141 give 

where Q. := WUW-^sj^rpy Similarly one has ||</'(0)||^cx)(jj) ^ ilo := l|f^(0)||7^s ^ • 
Observe that Assumption 12.11 implies that the conditions in Assump- 
tion are satisfied. 

With these preliminary remarks in hands, Theorem 14.51 vields 

with 

^{T) ■= ||(/i,Hf(^)'/2,Hf(^))|Il1{0,T;H1^) + 1 1 •^3,Hf(^) 1 1 Li (0,T;Hi) ' 
^'(T) := ||(/i,HF'/2,Hf)|Li(o,T;//1J + IU3,HF|Li(o,T;Hi)' 

One has ||t/||-^0(-j^) = — Jeu)U\\-^s ^^j'-^ and similarly ||t/(0)||^o^ = 
11(1 — Jei/)[/(0)||^s ^ f^o • Therefore, in view of the elementary inequal- 
ity X + y ^ 2x6^ (for all x ^ 1 and y ^ 0), the proof of Proposition 15.11 
reduces to establishing that ^{T) + ^'{T) ^ ^/TC{^l) . 

The estimate ^{T) ^ ^/TC{^^) immediately follows from the preliminar- 
ies. Indeed, define R and R' as in (j5.22|) . By definition of u := ^/JT+H, one 
has (J, ^ ^/Jlv and k ^ ^/kv . Hence, it is easily verified that Lemma 15.111 
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and Lemma l5 . 141 imply 

(5.24) ||/r;^p(^)||^, +||/2:hf(^)Li +\\f3UU)\\^,^C{R){l + R'}. 

Since C(-) is non-decreasing, integrating and using the Cauchy-Schwarz 
estimate, there results 

L2(0,T) /• 

This in turn implies the desired estimate 5(T) ^ \/TC{VI) since, by defini- 
tion, one has « ll^lli,o°(o t) + II-^'IIl2{o t) • 

Let us prove that, similarly, ^'{T) ^ ^/TC{^) . To do that it is sufficient 
to prove that (|5.24j) holds true with f^'^piU) replaced by /j'hf- This in 
turn follows from direct estimates. Indeed, observe that 

for some C°° functions J-i vanishing at the origin. As already seen, one can 
give estimates for the coefficients J-i by combining the Sobolev embedding 
Theorem with the Moser estimate (|5.7|) . It is found that 

(5.25) ||(/(,hf,/2,hf)IUi +I|/3,hf|Ihi ^^(^)ill^V^llHi + IIVMV^r||^i }. 
We next use Q^u ^ A'^ , to obtain 

IIv^^^IIh; + llv^^'^^ll/fi ^ + a/a* 11^11/^1+2 ^ R'- 

Consequently, the left-hand side of ()5.25|) is controlled by C{R)R' . 

To conclude the proof it remains to show that ^"{T) ^ ^/TC{^). This 
is nothing new in that it follows from the estimates: 

r(T)< ||(Ti,T3)||^.(o,^^^.+1) 

^ VrC(||(0, V^Vf)||^^([o^2^]^^)){l + ^||V7;||^2(o_y.^.+i)} 

^ A/rC(||i?||^oo(o,T)){l + II^'IIl2(o,T)} 

^ Vfc{vt). 

We have proved Proposition 15.11 This completes the analysis of the high 
frequency regime. 

6. Low FREQUENCY REGIME 

This section is devoted to the proof of a priori estimates in the low fre- 
quency region, which is the most delicate part. 

Proposition 6.1. Given an integer s > 1 +d/2, there exists a continuous 
non- decreasing function C such that for all a = (e, fj,, k) £ A , all T £ [0,1] 
and all U = {p,v,9) £ C^{[0,T]; H°°{B)) satisfying fll^ . 

(6.1) WJ.MWhut) ^ C(17o)e(^+^)^(^\ 

where u := y/jr+H , J^o := ||f^(0)||-^s , := ||C/||^s (j^) (see Definition \^^[ ). 
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As alluded to previously, the nonlinear energy estimates cannot be ob- 
tained from the L? estimates by an elementary argument using differentia- 
tion of the equations with respect to spatial derivatives. For such problems 
a general strategy can be used. One first applies to the equations some op- 
erators based on {edt) ■ Next, one uses the special structure of the equations 
to estimate the spatial derivatives. 

This basic strategy has many roots, at least for hyperbolic problems (see, 
e.g., ^ ^1 ) • For our purposes, the key point is that the hyperbolic 

behavior prevails in the low frequency regime. Yet, in sharp contrast with 
the Euler equations ( /i = k = ) , the form of the equations 1)2. 1|) shows that 
the time derivative and the spatial derivatives have not the same weight. 

In particular, our analysis requires some preparation. We begin our dis- 
cussion in H6.1l bv establishing some estimates which allows us to commute 
Jeui^dt)"^ with the equations. The latter task is achieved in ^6.31 

With these preliminaries established, we can proceed to give an estimate 
for Ji;y{edtYU . The fast components Jev{'A.\vv,Vp) are estimated next by 
using an induction argument. To conclude, we give the estimates for the 
slow components 6 and curl v . 

For the sake of notational clarity, in this section we deliberately omit 
the terms Ti and 5X3 in the system (|2.1|) . Nothing is changed in the 
statements of the results, nor in their proofs. 

6.1. Non-isotropic estimates. The fact that the time derivative and the 
spatial derivatives have not the same weight is made precise by the following 
lemma (whose easy proof is left to the reader). 

Lemma 6.2. There is family | -Ba,a 

\a € A, a G N^, 1 ^ |a| ^ 2} uni- 
formly bounded in C°° (M^; M^^^) (where N = {d + 2f) such that for all 
a ^ A and all smooth solution (p,v,6) of (|2.H) . the function ^ defined by 

(6.2) ^ := {^p,^t^,Vtp) where ^ := (9, sp, sv), 
solves 

(6.3) edt^= Baj{^)dj^ + e{fi + K) dj{Ba,jk{^)dk^). 

We want to introduce an operator based on (edt) which has the weight of 
a spatial derivative. The previous result suggests introducing the following 
family of operators. 

Definition 6.3. For all e ^ , ^ and £ G N, define the operators 

where, recall from 1)3. that A^^f := (/ — (ei^)^A) . 

We will need the following technical ingredient. 
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Lemma 6.4. Given F G C°°(R") satisfying F(0) = and (Tq > d/2, there 
exists a function C(-) such that for all e G [0, 1] , all v ^ [0,2] , all T > , 
all vector-valued function U £ C°° {[0,T]; H°° (B)) and all N B m ^ , 

(6.4) ||Z-F([/)||^.„_^C( \KuU\\h^,-.)- 

Proof. To prove this claim, observe that {£dt)'^F{U) is a sum of terms of 
the form 

f{U){edtY'u^---{edtY^Up, 

with p and X]i<i<p^« = m. In this formula, / is a C°° function and 
ui, . . . ,Up denote coefficients of U . 

In order to estimate these terms, we use the following result (whose proof 
follows from Proposition \^.2\ bv induction): 

Let a = (ao, , . . . , Op) G N^+^ be such that Yl^=o ~- 1"^! ^ "^o , then 

p P 

\-HY\v-\\ < TT II A~"'i/-ll 

i=0 i=0 

where the implicit constant depends only on d,ao and a . 

Set f ■= f — /(O) . We apply the previous result with oq = 0, ai = ii 
(i ^ 1), Vo = f{U) and Vi = {edtY'm (i ^ 1). This yields 



\A;:r{fmedtY^m...iedtY-Up)\ 



^ (l + ||/(^)||//<^o) l|-^e,V'"l||//<To-<i ■ ■ ■ ll-^efj^^pll/f'^o-^p- 

Since = /, the estimate (jSIHI) implies \\J{U)\\hoo ^ C(||Z° . 
Which completes the proof. □ 

Now we are in position to prove that Zl has the weight of a spatial 
derivative. The following result states that, for m G N, Z^^'^ satisfies the 
same estimates as A™F(^') does (where i*" is a given function). 

Proposition 6.5. Let s > 1 + d/2 he an integer. There exists a func- 
tion C(-) such that for all a = {e, fi, n) £ A, all T > and all smooth 
solution ip,v,e) £ C°°{[0,T]; H°° (B)) of (EH) , if £ [f^ + K,2] then the 
function ^ defined by (|6.2)) satisfies 



(6.6) IK'^Hm-^ ^ c{\mHs-r) iMhs 
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Proof. We prove by induction on m G {0, . . . , s} that 

(6-7) IK.^Ls-^-i ^^di^ii/^-O' 



i iTi 1 1 ^ /" 1 1 iTi 1 1 \ 1 1 iTi I , 



(6.8) \KuHHs~^<c{\mHs-^)m 

Note that these results are obvious with m = . 

Assume the results H6.7() - H6.7() at order m < s. By definition, one has 
= Ajj/Z™^ (edt) . It thus follows from Q that 

Since ^ ;U + k one has + K)djA~^ < I (see ()3.4() ). Similarly one 
has K'^^ < / . The proof thus reduces to estimating terms having the form 
Z^j^ . More precisely, it is sufficient to prove that, for all integer 

m < s and for all smooth function B , one has 

(6.9) ||Z- {B{^)d,^)\\^^_^_, ^ Cdl^'ll^.-,)' 

(6.10) ||Z- {B{^)d,^)\\^^_„^_, ^ C(ll^ll^.-x) ll^ll^. , 

where C(-) depends only on a finite number of semi-norms of B in C°° . 

Firstly, writing B{^)dj^ as djF{'^) for some C°° function F such that 
F(0) = 0, and using Lemma 16.41 with do = s — 1 , we find 

||Z™ (i3(^')5,vI/)||^„,<||Z™F(vI/)||^,_,_^c( Y 

As a consequence, the estimate (|6.9|) follows from the induction hypothe- 
sis (|6.7|) . Moving to the proof of (|6.1())) . we begin with the Leibniz rule 

Let N 9 £ ^ m. Since s > 1 + d/2 and m ^ s — 1 , Proposition \'A.2\ applies 
with do = s — 1 , o"! = mi = m — £ and (T2 = m2 = £ . It yields 

||Z™ (5(^^)9,^) < 

Using Lemma 16.41 to estimate the first term in the summand, we get 

Using the induction hypotheses (|6.7j) - (|6.7j) . we prove (|6.1flj) . □ 

We next prove a commutator estimate with gain of a factor e . This 
technical ingredient is an analogue of the commutator estimate (|5.2j) [with 
Q = h] we used in the analysis of the high frequency regime. 
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Lemma 6.6. Given s > 1 + d/2 , there exists a constant K such that for 
all £ G [0, 1] , all V ^ [0, 2] , all T > , a// m E N such that 1 ^ m ^ s and 
all f,u£ C°^([0,r];//°°(D)), 

(6.11) ||[/,J,,(eatr]n||^.-™+i 

m—l m—l 

Proof. The commutator [/, Ji;i^{edt)"^]u is expanded to 

(6.12) [/, JeuKedtru + J,,[f, iedtr]u. 

The proof of ()6.11|) is based on the tools we developed in Section |21 
a) We first claim that 

(6.13) II [/, J,,]{edtrn\\^s.^+, < ev ||/||^. ||Z- n||^_ . 

Since s > l+(i/2 and 1 ^ m ^ s, Proposition l3. 1 ll applies with (/i, m, (Tq, cr) ^ 
replaced by [ev^ 1, s, s — m) . It yields 

II [/, Jsu]{edtr4^^_^ < ev 11/11^. ||A-(2— _ 

Since m ^ s, there holds A^J^^ ^ A~J" . Hence, we have the first half 
of (|6.13() . namely: 

(6.14) II [/, Jeu]{edtru\\^,_^ < ev ||/||^. ||Z- n||^,_ . 

The technique for obtaining the second half is similar. We first apply Propo- 
sition 13.111 with m = , to obtain 

(6.15) II [/, J..](e5i)'"n||^_+, < 11/11^. ||Z-n||^_, 

and we next multiply (|6.15|) by ev . 

b) Moving to the second term in (|6.12|) . we claim that 

(6.16) \\j,,{f,{edtT\ 



U rjS — m + 1 

<eY,\\zl,dtf\\^^_,_,Y,\\Ze.^^^\ 



m-l , ^ 

(6.17) [f,{edtr]u = eY,[^lA{{edtfd,f){edtr-'^ 



■'e,u "'\\fJs-'^-t • 

Starting from the Leibniz rule, we get 

m—l 

u. 

£=0 

Let ^ £ ^ m — 1. Since s > 1 + (i/2 and 1 ^ m ^ s, Proposition 13.21 
applies with do = s — 1 , o\= m\ = £ and (T2 = m2 = m — 1 — £ . It yields 

\\A;,^-^'{{edtYdtf){edtr-'-'u)\\^,_^ 

< \\A-,\edtYdtf\\Hs-^-, \\Aj"^-'-'Hedtr-'-'u\\j,.-,-,^-^-,, . 



■^Here, m refers to the index used in the statement of Lemma 13.41 
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By summing over all ^ ^ ^ m — 1 , we obtain from the definition of the 
operators Z'^^^ that || A^J^^-'^ [/, (e9t)''"]ii||^s_„ is estimated by the right 
side of (|6.16|) . To complete the proof of use the elementary estimate 



(6.18) 



(Tp^yl/ rrS — m+l 



< 

J-[s-m ^ 



I A: 



I JJs — m 



□ 



6.2. Notations. Let us pause here to set a few notations we use continually 
in the sequel, as well as to make some running conventions. 

Notations. 

Prom now on, we consider a time < T ^ 1 , a fixed triple of parameter 
a = (e,/^,At) E A and a smooth solution U = {p,v,e) e C°°([0, T]; iJ°°(D)) 
of the system ()2.1|) . The notations (p, ip and ^ are shorthand notations 
for 

(6.19) (!>:= {d,ep), tp := {9, ep, ev) and ^' := (V^, (9*^, VV'). 
Here and below, s always denotes an integer s > 1 + d/2 . We set 

1/ := ^/J^+~K, ^0 ■= \\U {0)11-^8 and := ||C/||^s(-y) , 

where the norms are defined in Definition 12.41 As in (|5.22|) . set 

R ■= \\{P^v)\\j^s+i + \\9\\jjs+i , 

(6.20) , 

R' := + K \\{divv,Vp)\\jjs + a/^u || Vv||^s+i + a/k || V6'||^s+i . 

With these notations, one has Q « ||-R|Iloo(o t) + II-^'IIl2(o T) • 
Conventions. 

To say that a smooth non-decreasing function C: [0, +00) [l,-|-oo) is 
generic is to say that C(-) is independent of T, a and {p,v,9). Given a 
generic function C(-) , we denote by C the positive constant 

(6.21) C:=C(Oo)e(^+")^(^). 

The factor e is of no consequence but makes some arguments work more 
smoothly. We will often use the simple observation that for all generic 
function C(-) , one has C(0o) + iVT + e)C{n) ^ C. 

To clarify matters, with these conventions. Proposition 16. II is formulated 
concisely in the following way: there exists a generic function C(-) such 
that 1 1 <7^^ 1 1 ^'p^ ^ C . 

6.3. Localization in the low frequency region. 
Definition 6.7. Given m € N, set 

(6.22) ^™ := Je.iedtr. 
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Notation 6.8. Denote by /{J = (/™lf' /^lF' fsLp) commutator of the 
equations and the operator A"" : 

/{?LF := [91 (0), '^""] + [9i{cl>)v, X^] • - ^ A--] e, 

f^LF ■■= [92{cl>),X^]dtv+[g2{^)v,X^]-Vv-fi[B2{cl>),X^]v, 

The aim of this paragraph is to estimate . To begin with, consider 
the case when m ^ 1 . 

Lemma 6.9. There exists a generic function C(-) such that for all m £N 
such that 1 ^ m ^ s , 

II/lfI|-H=-™ •= (/iXF'/^Lf) H-s-m+l + /^LF [7=-™ + ! ^ C{R){1 + R'}, 

where R and R' are as defined in (|6.20j) . 

Proof. To emphasize the role of the slow component ^ (as defined in (|6.19|) ). 
many bounds are given in terms of: 

^■.= \\^\\fjs-i and r := ll^llji^s . 

Directly from Lemma l5.14| one has ||9t^||j:^s-i ^ C{R) and ^ 
C{R){1 + R'} . Furthermore, directly from the definition of ^ (see (|6.19() ). 
one has IKV', VV')||/^s ^ || (p, v) + ^ R- Hence, 

(6.23) 7 = ll^'llj:^,^! ^ C(i?) and T = W^^Wj^s !^ C{R){1 + R'}. 

Note that Proposition 16.51 applies since the condition v ^ fi + k is fulfilled. 
Step 1: Estimate for /{"lp • a) We begin by proving 

(6.24) ||[5i(0),A'"]atp||^._+i ^C7(7)r. 

Starting from Lemma 16.61 we find that the left-hand side is bounded by 

m—1 m 

(6.25) {\\gimHs + Y. \K-^^t9i{4>)\\Hs-i-.}{Y.\KAedtP)\\Hs-.}- 

e=o e=o 
Since "if = {. . . , {edt)p, . . .) , the estimate 1)6. 6|) implies that the second factor 
in H6.25() is estimated by C{'y)T. Moving to the first factor, note that 
dtgi{4>) = F{4>^dt4') for some function F such that F(0) =0. As a 
consequence. Lemma 16.41 implies 

m—1 m—1 

(6.26) Pludtgi{4>)\\Hs-,-, ^ c[Y. \\zlA<t^,M)\\Hs-.-)- 

Again, since ^ = {(j), . . . , dt(p, ■ ■ ■) , the estimate (|6.5|) implies that the right- 
hand side in the previous estimate is bounded by C{'y). It remains to 
estimate ||5i(</') ||//s . To do so write 

\\~9imHs ^ C{Mh.) ^ C{\\{^l^,Vn\Hs-^) ^ C{^)- 
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b) We next prove that 

(6.27) • Vp||^._+i ^ C{j){\\Vp\\hs-^+T}. 

Lemma 16.61 implies that the left-hand side is bounded by 

m—l 

{\\£9ii<P)v\\Hs+ Yl \\^e,udt{egi{(l>)v)\\^,_,_,^ 



(6.28) 



£=0 



1=1 

By definition ip = (9,sp,£v). Hence, one can rewrite £gi{(f))v as G{'ip) for 
some C°° function G. Consequently, the first term in the above written 
product is the exact analogue of the first term in (|6.25j) with gi{(p) replaced 
by G{tp) . We thus obtain that this term is estimated by C{'y) . 
Next, using the very definitions of , write 

m m 

(6.29) IK- Mm-' = E A-iV(eatp)||^._, . 
Since A~j/ < / , this in turn implies 

m m—l 

Which, as in the previous step a), is estimated by G{'y)T . 
c) To complete the estimate of /{"lp , we establish 

(6.30) ^ ||[5i(0),A""]0||^,_„+i ^ C7(7,i2){||K0||^.+2 +r}. 
Parallel to (|5.17|) . we decompose the commutator Ke~^[Bi{cl)),X''^~\9 as 

(6.31) ^ [£:!(</>), A"-] A^ + ^[£:2(0,v</.), A"-] -ve, 

where £i{^) ■= xi{£p)(3{e) and £2{<P,y^) := Xi{ep)P' {0)V6 . 
Replacing dtP with Ke~^A9, the arguments given in a) yield 



- II [fi((/>), X^]Ae\\^s-^+, ^ C(7) Y W^iu («A0)|U-* • 
We split the sum in the previous right-hand side as 

m 

e=i 

The first term is bounded by ||k6'||^s+2 . As regards the second term, we 
proceed as in b). Namely, write 

m m m— 1 

Y'^lK.^ow^^^.e = Y (eKA-»5i0||^._, ^ PU^A\hi-' ' 

i=l 1=1 1=0 
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where we used ekK'^^/S. < (which stems from k ^ ^/K ^ v). Note that 
^' = (. . . , dt9, . . .) . Therefore, the estimate (|6.6|) imphes that the right-hand 
side of the previous inequahty is controlled by C(7)r. 

We now have to estimate the H^~^~^^ -norm of the second term in (|6.3H) . 
Since k ^ i' and k ^ 1 , Lemma l6 . 61 implies that this term is bounded by 



m—l 



E=0 £=1 

This in turn is bounded by C{'y){C{R) + T} since 

(6.32) ||i/£:2(</>,V(/>)||^. ^ C{R), u\\zl,^t£2{ct>,^(t>)\\H.-^-e ^ C'(7)r, 

(6.33) ||V0||^..i ^ 7, V£G[l,m], \\uZl,Ve\\^^_, ^ C{j). 

The first inequality in 1)6. 32() follows from (|5.7I) . The second inequality 
in (|6.32j) follows from (|6.1U|) since the term dt£2{4>,'^<P) can be written 
-F(^')V^ for some function F . The first inequality in (|6.33|) is obvi- 
ous. In order to prove the last one, as already seen, one can write 

<l|z,y^'||^,_,_(,„,, ^c(7). 

d) By combining ()6.24p with (|6.27|) and ()6.3U|) . we obtain 

so that the desired estimate ||/i"lfIIj7''-'"+^ ^ C'(i2){l -|- R'} follows from 
the observation (|6.23|) . 

Step 2: Estimate for /^lf ■ Note that one can obtain /™Lp from /{"lp by 
replacing p hy v and 9 by ev . Therefore, we are in the situation of the 
previous step and hence conclude that 

(6.34) \\[92{4'),^'^]dtv\\Hi-^+i ^ C(7)r, 

(6.35) \\[g2{(p)v,X"'] •V?;||^.-™+i ^ C(7){ ||V^;||^.-i + T}, 

(6.36) ^||[B2(0),^'"]w||^=-.n+i i^Ci^,R){\\efiv\\j^s+2 + T}. 
Step 3: Estimate for /™lf • Note that, for ah u G , 

(6.37) ll""!! H'<T ^ e^"^ ll^ll H-CT ■ 
Hence, 

\\[g:,{<p),x'^]dte\\jjs-m+i + k\\[b^ ,x'^]e\\ 

^ i II [53(0),^"^] + - II [B3{c^),x"']e\\^s-^+, . 

The first term in the right-hand side is estimated as in a) above [replacing 
dtp by e^dt9] and the second term has been estimated in c). 
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For technical reasons, the estimate for [g3{(l))v, X"^]-'V6 is somewhat more 
comphcated. We argue as in the proof of Lemma ESI Set / := g3{(j))v and 
u = VO . We begin by sphtting the commutator [/, A""] • u as 

P + Q:= Jsu[f, {edtr] -u+if, J,,] ■ {edtru. 

By combining (|6.14|) with (|6.15j) multiphed by , we find 

11^11^—+-^ ^ WIWh^ ^ C'(7,i2)r, 

where we used (|6.33() . 

Our next task is to show a similar estimate for P . To do so we decompose 
P into two parts: 

(6.38) P1 + P2 := {Jeuif, {edtr] -u-Je^iu- {edtr f)] + J,,{u- {edtr f). 

Let us prove that ||Pi||^s-m+i ^ C{'y) . In light of H6.37() . all we need to 
prove is that 

(6.39) ||Pi||^,_™+i ^ eC{-f). 
Repeat the proof of Lemma 16.61 to obtain 

m—l 

\\Pl\\^s-rn+i \\^e,u '^tf Wjjs-i-e \\Z^,u'^\\H=--i--e ■ 

The sum differs from the one that appears in Lemma [6.6l in that it is indexed 
by ^ ^ 1 instead of £ ^ . This fact allows us to write 

m— 1 

\\Pl\\^s-rn+l < e ^ \\Ziu{edt)f\\fjs-l-l \\Ze7u^ dtu\\^^_-^_^ . 
l=\ 

Let 1 ^ £ ^ m — 1 . Write {edt)f as F{"^) for some function F such 
that F{0) = 0. By combining Lemma 16.41 and the estimate (|6.5I) . we obtain 

\Ku {edt)f\\Hs-^-,^C{j). 

Moving to the estimate of Z^~^ dtu, use the very definitions of u = V9 and 
^ = {. . . , dt9, . . .) , to obtain thanks to 1)6. 5() 

We have proved (|6.39() . so to conclude it remains only to estimate the 
second term P2 in 1)6. 38(1 . This is accomplished using 

ll-fzll/^s-'n+i ^ IKII_f/g \\^e!iy ^ {edt)f\\^s^-m+i , 

as the reader can verify, yielding the bound ||P2 1 1 ^ C{'y,R)T . 

This completes the proof of Lemma 16.91 □ 

Note that in the case when m = , the previous method does not work 
as it stands, as the estimate 1)6. 13p is no longer correct. 
Now we give an estimate valid for all ^ m ^ s — 1 . 
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Lemma 6.10. There exists a generic function C(-) such that for a// m £ N 
such that ^ m ^ s — 1 , 

(6.40) ||/Lp||^s-i-m := II (/£lf) /™lf) ll/^s-™ + ||/3?lf||/^s-™ ^ C{R). 



The estimate (|6.4()|) is nothing new in that it can be deduced by following 
the proof of Lemma 16.91 One only has to use the following analogue of the 
calculus inequality in Lemma 16.61 

\\[f,jsu{edtr]4m-- ^ 

m— 1 m—1 

In the case when m = the sum X^q ^ is interpreted as . The index s—2—£ 
in the second term of the first set of parentheses is not a typographical error. 
It is of use to us for the estimate of dt£2{4'i'^4') where £2 is as in 1)6. 32(1 . 

6.4. The fast components. We give here the estimates for the fast com- 
ponents div V and Vp . We make use of the notations introduced in ^6.21 

Notation 6.11. For all integer m ^ s, set Um '■= X^U where recall 



m 
•t) ■ 



As a preliminary step towards the estimate of {divv, Vp) we estimate the 
W°(T)-norm of U^- 

Lemma 6.12. For all integer m ^ s , there is a generic function C such 
that ||C^m|l7^o(j') ^ C, where C is as defined in (|6.2H) . 

Having estimated the commutators f^ , this result can be deduced by 
following the end of proof of Proposition 15.11 (see ^5.3|) . We therefore only 
indicate the points at which the argument is slightly different. 

Proof. It readily follows from Notation 16.81 that {p,v,6) := {pm,Vm,Gm) 
satisfies the linearized system (|4.1|) with 

/i := f^LF + /(,LF where /(,lf := -^Vxi(ep) • (/5(e)V0„) , 

/2 := f^LF + /2,LF where /£lf := f^X2{ep){2DvmVC{e) + divz;^Vr/(e)}, 

/3 := /s^LF + /3,LF where ^lf := '^X3(ep)V/3(^) • V0^. 

recalling that we deliberately omit the terms T and eT in the system 1)2. 1|) . 
Applying Theorem 14.31 we get 

WUmWno^T) ^ c(Oo)e^^^''^ l|f^m(o)||„o_^ + c(i7)5(r) + c{n)d'{T), 

with 

:= ||(/!Xf'/2?Lf)|Li(o,T;//1J + IU3?LfILi(0,T;//1)' 

^'{T) := ||(/i,LF>/2,Lf)|Li(o,T;//iJ + IU3,Lf|Li(o,T;//1)- 
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The proof thus reduces to estabhshing that ||C^m(0)||7Y0 ^ C*(f^o) 

5(r) + d'iT) ^ VTC{n) . To fix matters, we concentrate on the hardest 
case when m = s . Note that the conclusion of Lemma 16.91 is the exact 
analogue of Hence, one has ^{T) ^ \/rC(0) . As in all that 

has to be done in order to prove "S'iT) ^ ^/TC{VI) is to check -y/KH V^s||_ffi + 
V^II^^sIIhi^ ^ + i?'} . To do so, using the definition of 9s = 

j£y{edtY9 , we first rewrite VOs as eJ£y{edtY~^Vdt9 . Next, by combining 
the estimate ||eJe;yti|| rri ^ ||</e;/?^|| h-i ^ ||<^e;/^i|| ra with the inequality JH ^ 
V and the definition of ^' = (. . . , 5^^, . . .) , we obtain 

Analogous computations lead to 

^fJ^l^VsWni^ < \\vs\\H^ < \\ZI-'{edt)v\\^, ^ \\Z|-'^^,. 

Hence, the desired bound follows from 1)6.6^ and H6.23() . 

The technique for estimating the initial data is similar. Indeed, one 
has IIC/sIIt^o^ < ^fjl^j 1 as the reader can verify, yielding the bound 



||C/s(0)||^o_^ ^C(i?(0)) = C(Oo)- □ 

We now come to the main estimates of this part. 
Notation 6.13. Define ||ti||;c<'(T) ll^llL°°(o,r;_H"'^-i) + ^ II^IIl2(o,T;_H"'^) • 



Lemma 6.14. Let U := {p,v,9) be a solution of the system: 
(6.41) 



92{(p)dtv + ^Vp - fiB2{(l))v = f2, 



[ g^imO + divv - kB3{<P)9 = h. 

If the Fourier transform of U is supported in the ball {|^| ^ l/eu} , then 
there exists a generic function C{-) such that for all a S [l,s] , 

\\P\\jCZ+\T) + lldiv^^ll^.(T) 

C \\{edt)p\\,c.^T) + C \\{edt) dwv\\^.-i,^s 

(6.42) ^ ^ 

+ C" IIpIIl°°(o,T;L2) + C\\9{0)\\jj^+i +eC{n) ||Atw||^^+i(2.) 

+ eC(l^)||(/i,/2)||^.(^)+i^C||/3|| 

Remark 6.15. This result is an elaboration of Lemma 2.10 in pp. 

Proof, a) For further references, we first prove three estimates. 

1) Let (To > d/2 and a £ [0, ctq] . There exists a constant K such that for 

aU ^ and for ah (ui, U2) G H"" x (D) , 

(6.43) ||^il'"2||_H-'^ ^ ll^lll/f'^O ll^2||/i'<T . 
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To prove this result, we use a standard Moser's estimate. Since ctq > d/2, 
the product maps continuously x to for all r E [— cro,cro]- 

Hence, one has 

IKl^i2||//<T := ||?il^i2|lH<T-i + Q \\uiU2\\Ha 

^ IKi|Ih-0 (II'"2||//<t-i + q\\u2\\ho) = ||^i|Ih-o \\u2\\h^^ ■ 

2) Let a £ [0,s] and q G [1,+cxd] . Given a C°° function F , there exists a 
generic function C such that for all g ^ , and for all u G i^(0, T; i?^) , 

(6-44) \\FWu\\l^(0,T;H-) ^ C \\u\\l,^o^T;H-) > 

where ip is as defined in H6.19() . 

In light of (|6.43p . to prove this estimate we need only show that 

(6.45) \\Fm\L^io,T-Hs)^C. 

This will be established (independently) in (|().66)) below. 
2') Let us infer from (|6.44|) that, for all a G [0, s] , 

(6.46) ll^(V')^i|L-(r) ^ C , 
To see this write 

+ \\Fi^P)u\\ 

^ C'{||n||^oo(o,T;_H'<^-i) + ll^llL2(o,T;_ff^)} ^ ^ \\'^\\kZ{T) ■ 

3) Let cr G [0, s] , / G H'^-'^{B) and 7: D ^ (0, +00) be a function bounded 
from below by a positive constant. Furthermore, suppose 7 := 7—7 G -ff*(D) 
for some constant 7. We claim that if G H^(B) satisfies div(7Vn) = /, 
then there exists a constant K = K(d, s) such that 

(6.47) ||^i|li^<.+i ^ K h'^ll^oo {II/IIh— 1 + WiWh- + II'"IIl2 • 

Commuting the equation div(7Vii) = / with A°" and using the commutator 
estimate (|3.10|) to bound div([7, A'^]\7u) , we see that the proof of (|6.47|) can 
be reduced to the special case (7 = 0. On the other hand, the case o" = is 
immediate by usual integration by parts and duality arguments. 

b) Hereafter, RHS denotes the right-hand side of (|6.42l) . To prove that 

(6.48) ||divw||^^(y) ^ RHS, 
we begin by showing that 

(6.49) ||diviy||^j(^) ^ C{||(e5j)p||^.(2.) + ||K0^||^.+2(y) + e 
Rewrite div-u as 

-gi{ct)){edt)p + k£i{<I))M + k£2{<P,V<P)V9 + efi, 
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where £i and £2 ai'e as in (|6.31jl . Using (|6.46|) . one has 

Il5i(0)(e9t)p|lyc-(T) ^ ^ll(e50p|l/cj(T) ' 

Therefore, to infer H6.49() it remains only to estimate the /C^(r)-norm 
of k£^2(</', V(/>)V0 . Using the very definition of <S2(</', V(/)) and the esti- 
mate (|6.45|) . one has 

||'f2((/>, V(/>)||^oo(o,r;iT=-l) ^ C". 

Since s — 1 > d/2 , the estimate 1)6. 43p apphes with cio := s — 1 . It yields 
(6.50) 

||'?2(0, V(?!»)V0||^oo(o,T;/i'<^-i) ^ ||V0||^oo(o,T;/i'<^-i) ^ ^ II ^ (0,T;//'^) ' 

Moreover, one can easily verify that 

ly ||^2('^, V(/>)V^||^2(o,T;H-) 

^ lk'?2(0, V0)||^2(o,T;//=) \\'^^\\l°°{0,T;H'') 

^ a/T ||<?2(0, V(/))||^oo(o_T;Hg) II^IIl°=(0,T;_H''^+1) 

L°°{0,T;/i'J+^)) ll^llL°°(0,T;/i'<'+i) 
^ \/rC(r2) ll^ll Looj-Q^T./^CT+i) ^ C \\6\\ ioo(o,T;H<'+i) • 

From this together with (|6.5()j) . we conclude that 

11^2(0, V0)V6'||^a(T) ^ II^IIl°°(0,T;//<'+1) • 

This completes the proof of (|6.49|) since ||^||ioo(o t /Z'^+i) ^ II^IIa:'^+^(t) • 
So to prove H6.48|) it remains only to show that 

(6.51) ll«^^ll/cj(T) ^ R'HS. 

To see this, solve the first equation in (|6.41j) for div^; and substitute the 
result in the third equation, to obtain 

(6.52) 93(^0- K(i?3(0) - Bi{^))e = h + /s, 
with /g := gi{(t)){edt)p - efi . Since 

53((/)) - i?i(</.) := {x^iep) - xi(ep))div(/3(0)V-), 

the assumption Xiip) < Xsip) (which is (A3) in Assumption 12. ip implies 
that the equation (|6.52j) is parabolic. Hence, one can use the following 
estimate (see (|A.2() below): 

II'^^IIa:;J+^(T) ^ II V^i^^|lLoo(o,r;_ff<^+l) + IIV^^^^IlL2(o,T;_f/<^+2) 

(6.53) " _ ~ ' ' 

^C{||^Z.^(0)||^.+i+I.||/3 + /^|| 

From which we easily infer the desired bound (|6.51j) . 
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c) Estimate for Vp. Note that, by combining (|6.46|) with (|6.47jl . one has 

(6.54) W'^PWjcziT) ^ C'||div(52(0)~^Vp)||^^-i(y) + C||p||^^(o_7,.^2) . 
Starting from Vp = —g2{(t)){sdt)v + efiB2{(j))'v + e/2 , it is found that 

div(g2(0)""'^Vpi) = —{edt)divv + efi!Fi{(l))V^divv 

(6.55) „ „ , 

for some C°° functions J^i , J-'2 and J-^ vanishing at the origin. 

The most direct estimates show that the /CJ^"^(T) -norms of the first and 
the last three terms in the right-hand side of ()6.55|) are bounded by RHS. 
So we need only concentrate on the second term. Since ^ o" — 1 ^ s , the 
estimate (|6.46|l implies 

||£/iJP"i(0)V^ div 

The following observation is important since this where the spectral local- 
ization of U enters. Since v = Je^/sv and since e^V^ J^j^/s < , one has 
^ll/c°'^^(T) ^ ll'ii^ ^IIa:'^(t) ■ piece of information already de- 
termined in 1)6. 48() implies that this in turn is ^ RHS . Hence, the right-hand 
side of (|63i|) is ^ RHS . The proof of (|02|) is complete. □ 

Recall that the upshot of the two previous lemma is to estimate the fast 
components div Jgjyf and VJsvP ■ 

Lemma 6.16. There exists a generic function C(-) such that 

(6.56) \\JeuP\\L°°{0,T;H'') + ^ \\'^£'^P\\l'^{(},T;H^+^) ^ 

(6.57) ||div JeuV\\i^oo^Q^rp.jjs-l) + I' ||div -^£^^^11^2(0^7^;//=) ^ ^■ 

Proof. To avoid repetitions, we just give the scheme of the analysis. Set 

•= l|Pm||;cs-™+l(7i) + l|div Vm||;^s-m^j,-j . 

Setting U := Um and / := /{J — G{(p)v ■ VUm, where G((/>) is as defined 
in (|5.2.S|) . we are in the situation of the previous lemma and hence the 
estimate (|6.42() easily implies Xm ^ CXm+i + Ym where 

Ym ^ C{ ||f^m|l7-^0(T) + II (0) ||-^s-™ } 

+ {e + VT)C{0,){\\f[^\\^^^Q^rp.^s-l-m^ + ||C/'m||-^s-m^y-|}. 

Gathering the results of the previous lemma, one obtains 

y„ ^ c{c + c{no)} + {e + VT)c{n){c{n) + c{n)} ^ c. 

Hence, we end up with Xm ^ CXm+i + C . By an elementary induction, 
we therefore obtain Xq ^ CXg + C . Finally, noting that Xg < 1 1 1 1 7^0 (■p) 
and using Lemma 16.121 leads to the desired bound Xq ^ C . □ 
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6.5. The slow components. Having proved the estimates for the fast com- 
ponents, we now prove the estimates for the slow components 9 and curl v . 
To prove both estimates we use Assumption 12.21 This furnishes us with a 
C°° function S = S{^, p) such that 5(0, 0) = and 

(6.58) dS{{}, p) = gsi^, p) 6^ - g,{^, p) dp, 

(6.59) (■!?, p) {S{'&, p), p) is a change of variables. 

[Recah that p) G is the place holder of (6*, ep) .] 

Notation 6.17. Set a := S{9,ep) . 

One reason it is interesting to introduce the coordinate S = 5(i?, p) is 
that fj is well transported by the flow (see also Eemarks 16.201 and I6.22j) . 

Lemma 6.18. Given F € C°°(M) satisfying F{0) = 0, there exists a 
generic function C(-) such that 

(6.60) IK^t + t;-V)F(a) 11^.(0^^^^.) ^C(O), 

(6.61) ||F(a) 11^^(0^^.^.) < a 

Proof. Firstly, we form an evolution equation for a . Directly from the 
identity (|6.58|) . one has dt^x(^ = 9?,{4>)dt,x9 — £gi{4>)dt,xP ■ By combining the 
first and the last equations in (|2.1|) with this identity, we compute 

(6.62) dta + vVa- K{xz{ep) - Xii^P)) div(/3(0)Ve) = 0. 
Therefore, one has 

dta + vVa = KGi{()),Ve) + kG2(0)A0, 

for some C°° functions Gi and G2 , with Gi (0) = . Hence, the Moser 
estimate 1)5. 8|l implies at once 

\\dta + vVa\\Hs ^ C(||0||^.){l + l|V0||j^=+i}. 

Since ^/K ^ v , one has 

K \\V6\\^s+2 ^ \/k ||V6'||^s+i , 

so that \\dt(7 + V ■ Va\\jjs ^ C{R){1 + R'} where R and R' are as defined 
in (|6.20j) . Since Q, ||-R|Iloo(o t) + II-^'IIl2(o T) ' have 

(6.63) ||(5i + z;-V)a||^2(o,T;H=) ^C(l^). 

The chain rule and the rule of product in Sobolev spaces imply that the 
left-hand side of ()6.6U|) is estimated by 

(6.64) K{1 + ||F'((7)||^^(o,r^^.)) II (at + V ■ V)a\\ 

where F'{z) = F'{z) — F'{0) where F' denotes the differential of F . The 
first term in (|6.64|) is bounded by G{\\a\\j^cx,(^orp-H^)) ^ C{^1) (here is where 
we use the hypothesis that 5(0,0) =0). Hence, (|6.6fl|l follows from (|6.63|1 . 
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Moving to the proof of (|6.61|) , we first recall an usual estimate in Sobolev 
spaces for hyperbolic equations^: For all functions u,x ^ C^i[0, T]; H^{D)) , 
with s > 1 + d/2 , there holds 

(6.65) ||u||^..(o^r;H^) ^ h(0)|li^= + ^ e^^"*)^ \\dtu + x ■ Vn||^. dt, 



where X := K WxWh'^ some constant K depending only on (s, d) . 

The Cauchy-Schwarz inequality readily implies that the right-hand side is 
estimated by 

Applying this bound with (u, x) = (-^(<7), w) , it is found that the L°°{0, T; H'^) 
norm of F{a) is estimated by 

e™{C(17o) + Vt \\dtF{a) + v ■ VF(a)||^2(o,T;/f=)}- 
We complete the proof of ()6.61() by using ()6.6U() . □ 

Corollary 6.19. Given F G C°°(M?+'^) satisfying F(0) = 0, there exists a 
generic function C(-) such that 

(6.66) \\Fm\L^io,T;Hs)^C. 

Proof. The property (|6.59j) implies that there exists a C°° function F* such 
that F{'&, p, v) = F*{S{i}, p), p,v) for all (i?, p, v) E M x M x R'^ . Moreover, 
the hypothesis S{0,0) = implies F* (0,0,0) =0. 
Decompose F{7p) := F{6,ep,ev) = F*{a,ep,ev) as 

F*{a,0,0) + {F*{a,ep,ev) - F*{a,0,0)}. 

The first term is estimated by way of (|6.61|) . To bound the second term, ob- 
serve that one can factor out {ep, ev) . Consequently, there exists a function 
C(-) depending only on F* such that 

\\F*{a,ep,ev) - F*(a,0)||^. ^ eC{\\{e ,p,v)\\h.) ^ eC{R) ^ e^^^l'^H^^co.^)). 
Which completes the proof. □ 
Remark 6.20. The estimate (tPT^ implies ||9tF(^/;) ||^2(o,T;H=) ^ C{n). 
As in it yields ||i^(V')llL°°(o,T;H;) ^ C{no) + VTC{n) ^ C. Yet, this 

is weaker than (|6.66|) [indeed Hq = H^~^^]. 

Now we really use the special feature of the low frenquency analysis. The 
following lemma states that Ajj/cr satisfies parabolic-type estimates. 

Lemma 6.21. There exists a generic function C(-) such that 



^See also the estimate IIA.211 with 77 = . 
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Proof. In light of the estimate ||c||^oo(o t-H") ^ ^ (^^^ (IG-bip ). it suffices to 
prove that a := uA^JVa satisfies 

(6.67) II'''IIl°=(0,T;_H's) + II'''IIl2(o,T;_H's+1) ^ ^■ 

Let us form a parabohc evolution equation for a := uAjjj'S/a . Writing the 
identity (|6.58|) in the form di? = ci(t?, p)dS + 02(1?, p) dp with ci := l/g^^ 
and C2 := 51/53 , yields 

(6.68) Ve = ci (0) Vo- + ec2 {(l))Vp. 
Inserting this expression for V0 into the equation ()6.62() . yields 

(6.69) dta + V ■ Vcj — Hik{(j))/S.a = kGs + KG4, 
where k := (^3 — Xi)Pci is a smooth positive function^ and 

G3 + G4 := G3((A, V</)) + eG4(0, V0)Ap, 

for some C°° functions G3 and G4 , with ^3(0) = 0. 
We compute 

(6.70) dta + vVa - Kk{<t))/\d = g, 
where the source term is given by 

G := -i^Aj^iVv ■ Vfj) + i^[v, A-^] ■ VVa 

+ KAjXVk{(l))Aa) + Ki/[A:(0), A-i]AVa 

+ i^kA-^VG-s + iyKA^^VG4. 

Let us give the scheme of the analysis. In light of the standard esti- 
mate (|A.2|) (given in the appendix below), to prove (|6.67|) we need only 
show that the source term Q can be split as Qi + ^/KQ2 with 

\\G^\\l^{0,T;H^) ^ and \\G2\\l^(o^T;H^-^) ^ C*. 

To see this, there are only two nontrivial points. Firstly, note that 

||z/A/KA^J-VG4||^2(o,T;/i'=-l) 

:= \\iy^Aj^V{eGi{(p, V(p)Ap}\\ l2^o^t;H^-i) 

< ||V^G4(0, V(/>)Ap||^2(o,T;H=-i) (since ei^VA^^ < I) 

< (1 + ||G4(0, V(/))||^oo(o,T;H=-i)) II\/^pIIl2(o,T;H=+i) (straightforward) 
<CC. (hv (IMl .ljraill 

Secondly, there is a little trick: set 1^^^ := v j i^^fji^ \fK) ^ 1 and decompose 
vSIv ■ V(J into two parts: 



^By hypotheses, one has Xi < X3 ) < Ci ~ l/gs and < /3 (see Assumption 12. H . 
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Using the usual Moser estimates as well as the weighted versions (|5.6jl - (|5.7|l . 
it is easily found that 

Vl^\\VvVa\\Hs<\\Vv\\Hs WVaW^s < (i? + i?')C(i?), 

\\VvVa\\Hs-, < \\v\\Hsh\\j,s<RC{R), 

so that 

^/JI WVv ■ Va\\j^if^Q^j,.j^s) + \\Vv ■ V(t||^2(o,t;_h'''-i) ^ VtC{Q) ^ C. 
With these estimates established, the proof easily follows. □ 

Remark 6.22. One interesting feature of the equation H6.69() for a is that 
it is coupled to the momentum equation only through the convective term. 
Indeed, for the purpose of proving estimates independent of n, we cannot 
see the term div v [in the equation for 6] as a source term. 

We denote by curlw the matrix with coefficients (curlw)jj := djVi — diVj . 
The basic idea of the forthcoming computations is to apply the curl operator 
to the equation for v so as to cancel the large term e~^Vp. Yet, this 
requires some preparation because the factor g2{9,ep) multiplying the time 
derivative of v admits large oscillations in 0(1) . To get around this, we 
follow the analysis of • Namely, as in the proof of Corollary 16.191 we 
decompose g2{0,ep) into two parts: the first which is well transported by 
the flow, the second which admits small oscillations of typical size 0(e) . 

In particular, we do not estimate curlu directly. Instead we estimate 
curl(7of ) where the coefficient 70 is defined as follows. 

Notation 6.23. By one can write 52(1?; p) = ^{S{'&, p), p) for some 

smooth positive function T. Set ro(t?, p) := r(5('!?, p), O) and 

(6.71) 70 :=ro(^,ep). 

Lemma 6.24. There exists a generic function C(-) such that 

(6.72) ||curl(7ou)||^oo(o,T;H'=-i) + \^\\cu^Kiov)\\l^o,T;H-) ^ C. 

Proof. We begin by computing the equation satisfied by w := 0^1(7011) . It 
follows from elementary calculus that there exists a C°° function Ti such 
that for all (p, i?) G 

(6.73) ro(i9, p)/52(^?, P) = 1 + pri(i9, p). 

We first insert the expression for g2{0,ep) given by ()6.73() into the equation 
for V , thereby obtaining 

loidtv + V ■ Vv) + e~^Vp — fiB2V = F\ 

with F'l" := —pjiVp + eiJ,p'yiB2V where 71 := Ti{d,sp). Consequently, the 
equation for v is equivalent to 

{dt + v- V)(7ow) + e'^Vp - HB2V = F^ + {dao + v ■ V-fQ)v. 
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Using the elementary identity curlV = 0, we find 

{dt + V ■ V) curl(7ou) — /i-B2 curlu = 

where 

:= curl — [curl, v] ■ Vijov) + fi [ curl, B2] v + curl(u(9f7o + f • V70)) . 

Next, write B2 curl = 7o "^70-82 curl = 7o ^-62 curl(7o-) + 79 ^ [70, B2 curl] to 
obtain that uj = curl(7oi;) satisfies 

(6.74) dtUJ + vVuj- fi%^B2UJ = F, 

with F := Ftt + fi%'^[^o, B2 cm\]v . 
To sum up, F = Ylii<i< with 

Fi := - curl(p7iVp), F4 := ^[curl, ^2]^, 

F2 := e/ucurl(p7ii?2^^), F5 := -[curl, • V(7oi'), 

F3 := curl(t;(9j7o + v ■ V70)) , Fg := ^7o"^[7o, B2 c\xi\]v. 

Estimate for F . As in the proof of Lemma 16.211 to prove 1)6. 72() it suffices 
to show that one can decompose F as /i + y/Jif2 with 

(6-75) ll/illLi(o,r;i^»-i) ^ and ||/2|Il2(o,t;//=-2) ^ C"- 

To do so we decompose F as /i + with 

(6-76) ||/i|Il2(o,T;H=-i) ^ and ||/2|Iloo(o,t;H=-2) ^ C'(f^)- 

[Note that implies (IFTT^ since \/rC(0) ^ C.] 

Set /2 := e^p7i curl(S2f ) - ^%^v A 52V70 and /i := F - /2 . It 
follows from the very definition of that /2 satisfies the second estimate 
in 1)6. 76p . As regards f\ , the key point is the following: Starting from 

llFsll^,-! < [[wll^s ||9t7o + V ■ V7o||j;^s , 

we obtain ||F3||^2(o,T;iy»-i) ~ lbllL-'(o,T;H=) II^Ho + • V7o||i2(o,T;H=) • ^s a 
consequence, using the very definition of 70 (see 1)6. 71|) ). the estimate (|6.60|) 
implies that ||F3||^2(o,T;H=-i) ^ C{Vt). 

Let us estimate the other terms without repeated uses of the rules of 
product. To do so set P := {p,Vp) , V := {v,Vv) and 9 := {e,Ve). 
Direct computations show that one can write Fi + F2 + F4 + F5 + Fg — /2 
as 

(6.77) 5^(^5,.)"Ga,a(F, V, 9) (a G |a| ^ 1), 

for some C°° functions Ga,a such that the family {Ga,a | a G ^, |a| ^ 1} is 
bounded in , with in addition Go,a(0) = 0. The Moser estimate (|5.9|) 
implies 

\\{^d^rGa,a{P,v,e)\\jjs., ^ ||G,,a(p,y,9)||^. 

<c(||(F,y,9)||i..)||(F,y,9)||^. 
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Hence, directly from the definition of $7 , we arise at the desired bound 

Which completes the proof of the claim H6.76() □ 

6.6. End of the proof of Proposition 16. IL In order to prove Proposi- 
tion it remains to estimate Jeu(^ and JeuV ■ 

Lemma 6.25. Given F £ C°^(E^) such that F{0) = 0, there exists a 
generic function C(-) such that 

Remark 6.26. With = 9 this is the expected bound for JeuO . 

Proof. In light of (|6.61|) . it suffices to prove 

(6.78) l^\\JeuyF{(l))\\L^^^^T.jjs)+l^V^\\JeuV^F{(P)\\^2^0,T-,H^) ^ ^• 

By ()6.59|) one has 

(6.79) VF{(I)) = Fi {(l))Va + eF2 (0) Vp, 

for some C°° functions Fi and F2 ■ Applying Proposition 13.21 with ao = s, 
o"! = o"2 = mi = and m2 = 1 leads to 

||z.VA,-iF(</.)||^. < (1 + ||Fi(0)||^.) ||Vz.A7>||^. 

+ {l + \\F2{mH^)\Nei^KuP\\Hs, 

where F := F-F{0). Using i^\\-\\hs ^ WMhi^ Kv ^ ^ and (IfOmi . 
there results 

lk^^£!/"'"-^('?^)llL°°(0,T;//=) ^ ^{\\^eu'^\\L°°{Q,T;Ht+'^) + II^'IIl°°{0,T;_H"'=) } • 

As a consequence, from Lemma f6.21l and (|6.56|) . we get the first half of (|6.78j) . 

namely: ||i^-'ei/VF(^)||^oo(o,T;H'=) ^ ^ . 

The technique for obtaining the second half is similar. Differentiate ()6.79p , 
to obtain: didjF{(f)) := + J^u with 

Ti := Fi{(t))didja + eF2{(l))didjp, 
Tn := d^F^{^)dja + ed^F2{<^)djV. 

Exactly as above, one has 

Z^V^ II Jei.J^/||^2(o,T;//=+2) ~ C'l^V^ ll^£jy"^^llL2{0,T;Hs+2) + IIpIIl2(0,T;//=+1) } 

^ cc. 
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In order to estimate !Fii , note that J^jj := G{(l), V0) for some function 
G such that G{0) = . Hence, one has 

^Vk\\^ii\\l^{o,T;H^) ^ ^\\G{4','^(t>)\\L°°{Q,T,Hi) (straightforward) 

< Vtc(||<a||^^(o,t;hj+1)) (by (EZD) 

< VTC{n) ^ C (by definition). 
This completes the proof. □ 

The next estimate finishes the proof of Proposition 16.11 
Lemma 6.27. There exists a generic function C(-) such that 

Proof. We need only explain how to combine all the previous Lemma, 
a) To shorten notations, given a function / such that / — / G H'^ for some 
constant /, we denote by ||/||^<t the norm ||/ — /Hj^^^ • The semi- norm 
|M|/r<j is defined in the same way. 

We claim that for all smooth positive function 7 , there exists a constant 
K such that for all (e, ly) and for all vector field u G H'^{I}) , one has 

|Uei^^i||F7s ^ KM'^^^^\ \\diV J^iyUll + ||cUrl(7u)|| rrs-l 

(6.80) , 

+ \\Je„u\\^2 + e 

where M := 1 + ||7~"^||j|^s + IItII ^=+1 • 

To prove this result, we start from the following estimate: for all cr ^ 0, 

(6.81) < ||divn||j:^<^ + ||curln||^a + \\u\\l2 ■ 

As curln = 7"^ curl(7u) — 7~^V7 A u, by the usual rule of product in 
Sobolev spaces, we infer that for all a G [0, s — 1] , 

||curln||^„ < (1 + ||7~^llij.){||curl(7M)||^„ + ||7||^, 1^11^4 

<M||curl(7u)||^. +M2||n||^.. 

Thus, the estimate (|6.81() turns into 

ll'"ll_H"'^+i ^ ^^{l|divn||j|^a + ||curl(7ii)||^^ + ||n||j:^^}. 

By induction on N 9 o" ^ s — 1 , it yields 

(6.82) ||n||^. < M2«{||divn||^,_i + ||curl(7u)||j:^._i + ||n||^2}. 

Furthermore, by 1)3. 11(1 applied with (m, do, a) = (1, s + 1, s) , one has 

\\cuil{'yJeuu)\\jjs-i ^ II J£,.curl(7u) 11^,-1 + \\[y, J^uMlj^s 
< ||curl(7n)||j:^^ + ez^ ||7||^,+i ll'ullj:^, . 

From the estimate z^||7||£fs+i ^ Il7lliys+i ^ M and 1)6. 82p applied with u 
replaced by JeuU, the previous fact implies 1)6. 80() . 
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b) Let 7 := 70 (as defined in (|6.71|) ) and let M := 1 + II7 ^\\^s + ||7||^s+i . 
We list our bounds: 

by Lemma |6i25J : II-^IIl°o(o,t) 

by (|6.57() : ||div J£,^i;||^oo(o,T;H=-i) ^ C", 

by LemmainSH: l|curl(7ou)||^oo(o,T;/fs-i) ^ C, 

by Lemma 16.121 : ll-^ei^'W IIl°o(o,T;L2) ^ C, 

directly from the definitions : lk''^llL°c(o T /f") ^ ^ C- 

Hence, we deduce from (|6.8Up that || Jej^i'l 1^,00(0 t h") ^ ^ ■ This in turn 
implies WJeuvWioo^^Qrp.fjs+i^ ^ C , indeed recall that || Jejyw||j|^<T+i < ||'^£i/^^|lH<T • 

c) By (|6.81() . (|6.57p and Lemma [6. 121 to prove -v/a* ||</ei/V?;||^2(o t-h''+^) ^ ^ 
we need only check that ||curlj£;^f ||j^2(q j^.^i^s) ^ C. For this purpose, 
write 

\\curlJs„v\\j^s < (1 + ho'^\\Hs){\\cnTl{-foJeuV)\\jjs + ||7olli/s+i IbllnJ, 

and 

||curl(7oJei.t;)||^. < ||curl(7ov)||^. + holl^s+i \Mhs , 

where we used the commutator estimate ()3.11() with (m, do, a) = (0, s+1, s) . 
These two facts and the estimate ||7o||j;,oo(o r-jjs+i^ ^ ll7o||ioo(Q ^.^fs+i) ^ 

C (see Lemma l6. 25(1 implies 

^||curlJ^^u||^2(o,T;H=) ^ C'{\//^l|curl(7of)||i2(o,T;H'=) + \H L^0,T;H-)} ■ 

Hence, the expected bound directly follows from Lemma 16.241 and the 
estimate ||'y|lL2(-o^T;//s) ^ \\''^\\l°°{o,T]H^) . □ 

7. Uniform stability 

In this section we prove Theorem 12.71 Granted the uniform a priori 
bounds proved in the previous sections, the proof essentially reduces to 
establishing a local existence result for fixed a G A. 

7.1. Local existence result. For fixed a £ A, the fact that the Cauchy 
problem for (|l.lj) is well-posed is immediate provided that one chooses to 
work with the unknown {p,v,T) . Indeed, with this choice of dependent 
variables the system is a coupled hyperbolic/parabolic system in symmetric 
form, so that the general theory applies. We begin by proving that this 
remains valid for the system (|2.1j) . 

Introduce first a definition clarifying the structure of the systems we pro- 
duce. Consider a system of nonlinear equations 

(7.1) Ao{u)dtu+ Aa{u)d^u = F{u,Vu), 
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where each Aj^ (0 ^ ^ 2) is a n x n matrix smooth in its arguments 
and furthermore symmetric: Afj = A^^ . We suppose F is smooth in its 
arguments, with values in M" . 

Definition 7.1. The system (|7.1() is said to be of coupled hyperbolic/parabolic 
type provided that there exists (ni, 71,2) G , and a splitting of the unknowns 
u = {u^,u^) £ M"i X M"2 50 ifiat F{u,Vu) = (0, ^^(n, Vn))* and 

^°(^)=V Af{u))^ ^"(^)=V0 ^^2(^)J /or a// |a|= 2, 
where in addition 

V(u,0 G X S'^-\ Al\u) > 0, Al^u) > 0, - ^ Nr > 0. 

\a\=2 

(here M^^ denotes the sub-blocks of the matrix M which corresponds to the 
splitting W3u = (n\ u^) G M"i x ). 

The general theory [281 l4Uj for hyperbolic or parabolic problems applies 
for coupled hyperbolic/parabolic ones. One can prove the following result. 

Proposition 7.2. PI) Let s > 2 + d/2. Assume that (|7.1() is of coupled 
hyperbolic/parabolic type. Then for all initial data uq in H'^{]D) there ex- 
ists a positive T = T{\\uo\\jjs) such that the Cauchy problem for 1)7.11) has 
a unique classical solution ^{uq) G C^{[0,T]; H^{B))) such that ^'^{uq) G 
L'^{0,T; H^~^^{B))) where ^'^{uq) are the last n2 components of ^{uq) where 
n2 is as in Definition \7.1\ 

P2) The interval [0,T*) , with T* < +00, is a maximal interval of 
existence if and only if limsupi_>j'* ||^(^^o)(i) |lvKi'°=(D) ~ • 

Remark 7.3. i. If the system is symmetric hyperbolic, it is enough to 
assume s > 1 + d/2. 

ii. Let u G C°([0, T]; i?''(D)) be a classical solution of dZH). It follows from 
property P2) that, if m(0) G H'^{B) then u G C°°([0, T]; if°^(D)) . 

Fix a = {£,11, k) G a. We next show that an appropriate change 
of variables transforms the system ()2.H) into a system of coupled hyper- 
bolic/parabolic type. Let g denotes the function given by Assumption 12.21 
Since the mapping 3 p) i-^' g{'d, p)) G is a C°° diffeomorphism, 
it is equivalent to work with {p,v,6) or with {p,v,6) where p := q{9,£p) . 

Firstly, we form an evolution equation for p := g{6, ep) . By combining 
the first and the last equations in ()2.H) with the second identity in (|2.2|) 
written in the form Xs^Q = X39i dp — Xi93 di? , we get 

(7.2) X3idtP + v-Vp) + {X3 - Xi) divv = 0. 

The second identity in 1)2.2(1 written in the form 

dp = 71 d-!? + 72 dg := _)_ J_ dg^ 

X391 91 
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yields eVp = jiVO + 72 Vp. Summing up, one has 
' XsidtP + v • Vp) + (x3 - Xi)divv = 0, 

(7.3) I g2{dtv + vVv) + e-^jiV9 + e-^j2^p - pBiV = 0, 

, gsidtO + v-Ve) + divv- KB39 = 0, 

where all the coefficients are positive by Assumption 12.11 We also mention 
that all the coefficients are evaluated at {9,ep) = {6, P{9,p)) for some C°° 
function P . We next symmetrize the system ()7.3() . To do so multiply the 
first equation by (x3 — Xi)~^7i > the second by and the third one by 72 , 
to obtain 

' 5i{dtp + V ■ Vp) + 71 divv = 0, 

(7.4) < 52{dtv + V ■ Vv) + 71 Vp + 72V6' - e'^pB2V = 0, 

, '^3(5*6' + • V6') + 72 divu - K72B36I = 0, 

with 61 := (x 3 - Xi)"^7iX3 , h ■= g^ ff2 a nd ^3 := 7253 [compare with the 
system 1)4. ISp ]. Again, by Assumption [23 the coefficients 5i 's are positive. 

Lemma 7.4. For all (/i, k) G [0, 1]^ , the symmetric system (|7.4j) is of 
coupled hyperbolic/parabolic type. 

Proof. Set n = 2 + d . Ifp = K = set ni = n and v} = u. If p = and 
K > set ni = 1 + d and = (p, f ) . If p > and k = set ni = 2 and 
= {p,9) . If p > and k > set rii = 1 and v} = p . □ 

With this result in hands, we obtain that the Cauchy problem for (|7.4() is 
well-posed for all fixed a S >1 , so is the Cauchy problem for ()2.1() since 1)7. 4() 
has been deduced from ()2.1|) by using a C°° diffeomorphism. 

We are now prepared to prove Theorem 12.71 

7.2. Proof of Theorem EH Let s > 1 + d/2 and Mq > 0. The previous 
analysis implies that, for all a G ^, there exists Ta = Ta{MQ) such that for 
all initial data Uq in the ball B{Ti.'^Q, Mq) , the Cauchy problem for (|2.1|) 
has a unique classical solution Ua in TC'^{Ta) ■ Denotes by T* the maximal 
time of existence of such a classical solution. For all t < T* , set i^a{t) '■= 
||C/a||-^s • task is to show that there exists T > and M < +00 such 
that 

(7.5) Va G A, T* ^ T and Qa{T) ^ M. 

Up to replacing the initial data Uq by JsUq and letting 5 goes to zero, 
we can assume that Uq G H°° (B) . In light of Remark 17. 3( it means that 
we can assume Ua G C°°{[0,T*);H°°{I})) . By combining Propositions 15. ll 
and 16.11 we know that there exists a smooth non-decreasing function C(-) 
such that for all a £ A and all t < min{T*, 1} , there holds 



(7.6) 



na{t) ^ C(Jla(0))e(^*+^)^(^''(*». 
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Choose first Mi > C{Mq) and next Ti ^ 1 and ei ^ 1 such that 
(7.7) C7(Mo)e(^+^i)^(^^i) < Mi. 

Let t < min{T*,Ti} and a = {e,fi,K) G A be such that e ^ ei . By 
combining the inequahties H7.6() and ()7.7() with the hypothesis ila(O) < Mq , 
we infer that i^a{t) 7^ Mi . Besides, we can assume without restriction that 
Mq < Ml , so that ria(O) < Ml . Since the function f^a is continuous, we 
infer ^la{t) < Mi . Consequently, the continuation principle [which is P2) 
in Proposition I7.2j shows that T* > Ti. 

On the other hand the [omitted] proof of Proposition 17.21 implies that 
there exists T2 > and M2 < +00 such that for all a = (e, ^, k) G A with 
e ^ ei , one has T* ^ T2 and ^^(rs) ^ M2 . 

We have proved ((73|) with T := min{ri/2, T2} and M := max{Mi, M2} , 
which completes the proof of Theorem 12.71 

8. Decay of the local energy 

We will consider systems which include (|2.H) as a special case. The mo- 
tivation is to consider a structure general enough to include, say, the com- 
bustion equations. To do so we allow the limit constraint on the divergence 
of the velocity field to read divv = F{D^tp) where V is the slow variable 
(namely dtip = 0(1) ), m ^ is a given integer and 

More precisely, we consider system of the form: 

gi{i^)dtP +-^d\wv = ^Fi(Z)™V) + i"2(^™V', D'^v), 
^^■^^ { g2{'4>)dtv+\vp = F^{D^i,,D^p,D"^v), 

where {p, v, ip) is defined on [0, T] x M*^ with values in M x x , n ^ 
and m ^ are given integers. 

Assumption 8.1. The functions gi (i = 1,2) and Fj (j = 1,2,3,4) are 
smooth in their arguments, satisfying 

5i(y),52(y) >0, Fi{Y),F2(Y,Z)eR, Fs{Y,Z)eR'', F4(y,Z)GM", 

where y , Y and Z are the place holders of ip , D'^ijj and {D"^p, D'^v) . 
Moreover Fj{0) = for j = 1,2,3,4. 

Assuming that the slow variable decays sufficiently rapidly at spatial in- 
finity we prove that the penalized terms converge to . 

Proposition 8.2. Let T > , m ^ and s > 2m + d/2 . Assume that the 
functions {p^,v^,ip^) satisfy ()8.1|) and 

(8.2) \\{p',v',r)\\^s^T) ■■= sup ||(/,^;^)(t)||^. + ||^^(t)||^.+_i < +00. 

te[o,r] 
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Assume further that i/j'^ converges strongly in C^([0, T]; H^^^{W^)) , for some 
a > 1 + d/2, to a limit tp satisfying, for all {t,x) G [0, T] x , 

(8.3) \xp{t,x)\ !^ K\x\-'^-^ , \Vip{t,x)\ ^ K\x\-^-^ , 

for some positive constants K and 7 . Then, converges strongly to in 
L'^{0,T; H^^^{W^)) for all s' < s and dw v"^ — Fi{D"^ip'^) converges strongly 
to in L2(0, T; H[^^{R<^)) for all s' < s - I . 

Proof. Since (p^ , , -0"^ ) is uniformly bounded in /C®(T) (see H8.2|) 'l. using 
an interpolation argument, it is sufficient to prove that p^ and divu^ — 
Fi(D™'0^) converges strongly to in L'^{0,T; Lf^^{R'^)) . To prove this 
result we use the following Theorem proved in [HJ [although not explic- 
itly stated in this way]. To clarify matters, we mention that: (1) all the 
functions considered in the following statement are real-valued; (2) all the 
convergences considered in this proof are strong. 

Theorem 8.3 (Metivier & Schochet). Let T > and let be a bounded 
sequence in C°([0, T]; ^^(R'^)) gu^h that 

e^dt{a^dtu^) - div(6"Vn^) = c^ 

where the source term converges to in L'^{Q,T; Lp'iW^)) . Assume fur- 
ther that, for some cr > 1 + d/2, the coefficients {a^,b^) are uniformly 
bounded in C^{[0,T]; H" {R"^)) and converges in C^{[0,T]; HI^^^{M.'^)) to a 
limit (a, 6) satisfying the decay estimate 

\a{t,x)-a\ ^ Klxp^"^, |Va(t,x)| i^K\x\-^~^ , 

\bit, x)-b\ ^ K \x\-^~^ , \Vb{t, x)\ \x\-^~"' , 

for some given positive constants a,b,K and 7 . 

Then, the sequence converges to in L'^{0,T; Lf^^{W^)) . 

We can directly apply Theorem 18.31 to prove the first half of Proposi- 
tion |SI21 that is the convergence of p^ to in L'^{0,T; Lf^^{R^)) . Indeed, 
applying e'^dt to the first equation in 1)8. we compute 

(8.4) e^dtia'dtp') - div(6"Vp^) = c' := ef, 
with := gi{ip^) , b^ := \/g2{ip^) and 

(8.5) r := dtFi{Y') + edtF2{Y',Z')-d\^{b'F^{Y',Z')), 

where Y'- := D"^^ and := {D'^p^ , D'^v^) . 

The equations (|8.1|) imply that one can rewrite as 

where .7-" is a C°° function such that J-{0) = . Using the usual nonlinear 
estimate in Sobolev spaces ()5.8p . the hypothesis 1)8. 2|) implies that is 
uniformly bounded in C°{[0,T]; H'"^-^ (R'^)) . Consequently, = ef^ 
converges to in L^{0,T;L^{R'^)) and Theorem lO applies. 
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To prove the second half of Proposition 18.21 we begin by proving that 
:= {edt)p^ converges to in L^(0, T; L^q^(R'^)) . To do so apply {edt) on 
the equation ()8.4() . to obtain 

with 

f := edtf - edt{dta%edt)p') + dw{dtb'Vp'), 

where is given by H8.5() . Again one can verify that is a bounded 
sequence in C°([0, T]; L^(M'^)) , which proves the desired result. 
To complete the proof, observe that 

divT;^ - Fi(y^) = -g^{^l^^){e^t)f + eF^iZ'). 

Hence, the fact that divw^ - Fi{Y^) converges to in ^^(0, T; Lp^^(M'^)) 
follows from the previous step and the fact that gi{ip^) — 51(0) and F2{Z^) 
are uniformly bounded in C°([0, T]; . □ 

8.1. Proof of Theorem 12.91 . To simplify the presentation we concentrate 
on the hardest case when k is a fixed positive constant. 

We first prove the convergences for some sub-sequence of {p^,v^,6^) . 

The equation for 9 implies that dt9 is bounded in C°{[0,T]; H^-^iW^)) . 
Therefore after extracting a sub-sequence, we can assume that, for all s' < s , 

(8.6) 6' ^9 inC70([0,r];/7f;+i(IR'^)), 

where the limit 9 belongs to C°([0, T]; Hf^+^ {R'^)) n L°°(0, T; H'+\R'^)) . 

Since {p'^,v'^) is uniformly bounded in {[0,T]; [W^]) , after extract- 
ing further sub-sequence, we can also assume that 

(8.7) {p',v') {p,v) weakly * in L°^{0,T; H'(R'^)). 

Note that the system (|8.ip includes system (|2.1() as a special case where 
:= (6^,ep^) and m = 2. It follows from the very definition of and 

the norms IMIfrCT that 

e 

\\ip',v')\\H. + W\\H.+m-i ^ \\ip',v')\\Hs+e\\p'\\Hs+i + 

^ll(/,^^)ll^.+i + lir||^.+, 
<\\ify)\\Hit' + W\\m+^, 

where the implicit constant depends only on the fixed positive value of := 
Vl^ + i^- It follows that ||(pS?^S'0^)||/c'=m ~ m- Asa 

consequence the first hypothesis of Proposition 18.21 is satisfied. Therefore, 
in order to apply Proposition 18.21 to the equations (|2.H1 so as to obtain 
Theorem 12.91 it only remains to show that = (^,0) satisfies the crucial 
assumption (|8.3|1 . 
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Multiplying the first equation in (|2.1j) by e and passing to the weak limit 
(in the sense of distributions) shows that divv^ — div(/?(0^)V6'^) 
converges to in that sense. Therefore, the limit 6 satisfies 

(8.8) gs{e, o){dte + vve)- ^(^3(0) - xi(o)) diY{p{e)ve) = 0. 

The assumption xs > Xi [see (A3) in Assumption 12. Ij implies that (|8.8|) is 
parabolic. The desired spatial decay thus follows from estimates in weighted 
Sobolev spaces. Indeed, introduce 6s '■= {x)^9 for some given 6 > 2. Then 

dtOs - dw{k{e)Ves) = Vc^it, x)d^es (« G N*^, \a\ < 1), 

for some smooth positive function k and some coefficients Va such that 
||Kt||j7s ^ C'dK^i ^)ll_f/s) • Commuting A'' with the equation and applying 
the usual estimate, we obtain 

(8.9) \\Os\\L^io,T;Hs)^Ce^''\\es{0)\\Hs, 

for some constant C depending only on || (w, ^) ||^oo(o y-i^s) (see ()A.2|) with 
r] = 1). The convergence 1)8. 6() implies that 9{0) is the limit of the initial 
data 61^(0) . Using the hypothesis (x)'^6'o £ H%W^) , it follows from (jHSj) that 

Since s > l + (i/2, the Sobolev Theorem implies that H%R'^) ^ VF^'°°(M'^) . 
Hence, the assumption (|8.i-i|) is satisfied. 

A moment's thought shows that the convergences holds for the full se- 
quence. The proof of Theorem 12.91 is complete. 

Appendix A. Estimates for hyperbolic and parabolic systems 

Here we briefly recall some standard estimates for hyperbolic or parabolic 
systems we used throughout the paper. We state estimates for the solutions 
u = u{t, x) £ of systems having the form: 

(A.l) dtu + ^ A^dy -r]Y, Aady = h + ^^/s, 

|a|=l |a|=2 

where each = Aa{t,x) is a n x n matrix valued-function smooth in its 
arguments and furthermore symmetric: Ap = Aj^ ; also, the source terms 
fi = fi{tix) are smooth functions with values in . In particular, when 
r/ = , the system (|A.1|) is symmetric hyperbolic. We assume that for > 
the system is parabolic: for all triple (t, x, ^) with ^ 7^ , 

A4t,x)r >o. 

I"l=2 
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Lemma A.l. Let ^ cr ^ s G (l + d/2,+oo) . There exists a function C(-) 
such that for all T G [0,1], all u G C^([0, T]; F°°(D)) satisfying (TO]) and 
all T] G [0, 1] , 

II^IIl°°(0,T;//'^) + II'"'IIl2(0,T;_H''^+1) 
^"^■^^ ^ (lln(O)ll^. + ||/l||il(o,T:H'^) + ll/2lL2(o,T;H-i)) , 

where K = K{s,d) and R := '^^^^^2\\^a\\L°°{0T-H^) where A = A — A{0). 

Proof. Commuting the equation <\A.l\i with A'^ , using the commutator es- 
timate (jc{.lU|l and the Gronwah's lemma, the proof of <\A.2\i can be reduced 
to the special case o" = . The latter case relies upon usual integrations by 
parts and duality arguments. □ 
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